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THE STABILITY OF THE LAMINAE BOUNDARY LAYER IN A COMPRESSIBLE FLUID 

By Lester Lees 


SUMMARY 


The present paper is a continuation of a theoretical investi- 
gation of the stability of the laminar boundary layer in a com- 
pressible fluid. An approximate estimate for the minimum 
critical Reynolds number R e , or stability limit, is obtained 

CT min 


product of the mean density p* and mean vorticity across the 


in terms of the distribution of the kinematic viscosity and the 

, du* 
dy* 

boundary layer. With the help of this estimate for Re CTmin 
it ■ is shown that withdrawing heat from the fluid through the 
solid surface increases R e and stabilizes the flow, as 


compared with the flow over an insulated surface at the same 
Mach number. Conduction of heat to the fluid through the 
solid surface has exactly the opposite effect. The value of Rscr mim 
for the insulated surface decreases as the Mach number in- 
creases for the case of a uniform free-stream velocity. These 
general conclusions are supplemented by detailed calculations 
of the curves of wave number { inverse wave length) against 
Reynolds number for the neutral disturbances for 10 representa- 
tive cases of insulated and no-ninsulated surfaces. 

So far as laminar stability is concerned, an important differ- 
ence exists between the case of a subsonic and supersonic free- 
stream velocity outside the boundary layer. The neutral 
boundary-layer disturbances that are significant for laminar sta- 
bility die out exponentially with distance from the solid surface; 
therefore , the phase velocity c* of these disturbances is subsonic 
relative to the free-stream velocity u a * — or u 0 *—c*<i. aff, 


where aff is the local sonic velocity. When AA=M 0 <f 1 

{where M 0 is free-stream Mach number), it follows that 
0 ^ c* ^ c * m ax; on d any laminar boundary-layer flow is ultimately 
unstable at sufficiently high Reynolds numbers because of the 
destabilizing action of viscosity near the solid surface, as 
explained by Prandtlfor the incompressible fluid. When M 0 f> 1 , 


however, ==>1- 


1 


> 0 - If the quantity 


d ( —du* 

Jy*V' 


dy*/Jp=c 


is large e nough negatively, the rate at which energy passes from 
the disturbance to the mean flow, which is proportional to 
- du*'' 


Sd f—du*\ 

c U/*\ p 


j— * ) _ , can always be large enough to coun- 


terbalance the rate at which energy passes from the mean flow 
to the disturbance because of the destabilizing action of viscosity 
near the solid surface. In that case only damped disturbances 
exist and the laminar boundary layer is completely stable at all 


Reynolds numbers. This condition occurs when the rate at 
which heat is withdrawn from the fluid through the solid surface 
reaches or exceeds a critical value that depends only on the 
Mach number and the properties of the gas. Calculations 
show that for M<R>8 {approx.) the laminar boundary-layer flow 
for thermal equilibrium — where the heat conduction through the 
solid surface balances the heat radiated from the surface — is 
completely stable at all Reynolds numbers under free-flight 
conditions if the free-stream velocity is uniform. 

The results of the analysis of the stability of the laminar 
boundary layer must be applied with care to discussions of 
transition; however, withdrawing heat from the fluid through 
the solid surface, for example, not only increases Re C r mia but 
also decreases the initial rate of amplification of the sel f-excited 
disturbances, which is roughly proportional to 1 j-\/Rt crmU . 
Thus, the effect of the thermal conditions at the solid surface on 
the transition Reynolds number Re lr is similar to the effect on 
R)a- min ■ A comparison between this conclusion and experi- 
mental investigations of the effect of surface heating on transi- 
tion at low speeds shows that the results of the present paper 
give the proper direction of this effect. 

The extension of the results of the stability analysis to laminar 
boundary-layer gas flows with a pressure gradient in the 
direction of the free stream is discussed. 


INTRODUCTION 

By the theoretical studies of Heisenberg, Tollmien, 
Sehlichtmg, and Lin (references 1 to 5) and the careful 
experimental investigations of Liepmann (reference 6) and 
H. L. Dryden and his associates (reference 7), it has been 
definitely established that the flow in the laminar boundary 
layer of a viscous homogeneous incompressible fluid is un- 
stable above a certain characteristic critical Reynolds 
number. When the level of the disturbances in the free 
stream is low, as in most cases of technical interest, this 
inherent instability of the laminar motion at sufficiently 
high Reynolds numbers is responsible for the ultimate 
transition to turbulent flow in the boundary layer. The 
steady la min ar boundary-layer flow would always represent 
a possible solution of the steady equations of motion, but 
this steady flow is in a state of unstable dynamic equilibrium 
above the critical Reynolds number. Self-excited dis- 
turbances (Tollmien waves) appear in the flow, and these 
disturbances grow large enough eventually to destroy the 
laminar motion. 
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The question naturally arises as to how the phenomena 
of laminar instability and transition to turbulent flow are 
modified when the fluid velocities and temperature varia- 
tions in the boundary layer are large enough so that the 
compressibility and conductivity of the fluid can no longer 
be neglected. The present paper represents the second 
phase of a theoretical investigation of the stability of the 
laminar boundary-layer flow of a gas, in which the com- 
pressibility and heat conductivity of the gas as well as its 
viscosity are taken into account. The first part of this work 
was presented in reference 8. The objects of this investiga- 
tion are (1) to determine how the stability of the laminar 
boundary layer is affected by the free-stream Mach number 
and the thermal conditions at the solid boundary and (2) to 
obtain a better understanding of the physical basis for the 
instability of laminar gas flows. In this sense, the present 
study is an extension of the Toflmien-Schlichting analysis of 
the stability of the laminar flow of an incompressible fluid, 
but the investigation is also concerned with the general 
question of boundary-layer disturbances in a compressible 
fluid and their possible interactions with the main external 
flow. 


SYMBOLS 


With minor exceptions the symbols used in this paper are 
the same as those introduced in reference 8. Physical 
quantities are denoted by aD asterisk, or star; whereas the 
corresponding nondimensional quantities are unstarred. 
A bar over a quantity denotes mean value; a prime denotes 
a fluctuation; the subscript 0 denotes free-stream values 
at the “edge” of the boundary layer; the subscript 1 denotes 
values at the solid surface; and the subscript c denotes 
values at the inner “critical layer,” where the phase velocity 
of the disturbance equals the mean flow velocity. The 
free-stream values are the characteristic measures for all 
nondimensional quantities. The characteristic length meas- 
ure is the boundary-layer thickness 5, except where otherwise 
indicated. Note that in order to conform with standard 
notation, the symbol 5 for boundary-layer thickness is 
unstarred; whereas the symbols 5* and 8 are used for 
boundary-layer displacement thickness and boundary-layer 
momentum tliiekness, respectively. 


V* 

t* 


distance along surface 
distance normal to surface 
time 

component of velocity in redirection 


w—- 




r 

p * 
p* 
y* 


V * 


component of velocity in ^-direction 


stream function for mean flow 
density of gas 
pressure of gas 
temperature of gas 
laminar shear stress 


k* 

c$ 

Cp 

R* 

y 

c* 

X* 

S 

S* 


d 


ordinary coefficient of viscosity of gas 

kinematic viscosity of gas (pi*/p*) 

thermal conductivity of gas 

specific heat at constant volume 

specific heat at constant pressure 

gas constant per gram 

ratio of specific heats (c„/c s ); 1.405 for air 

complex phase velocity of boundary-layer 

disturbance 

wave length of boundary-layer disturbance 
boundary-layer thickness 
boundary-layer displacement thickness 

(J o (1— pw)dy*J 


boudary-layer momentum thickness 
pw(l—w)dy*J 


a * _ wave number of boundary-layer disturbance 

(2r/x*) 

2tt 

a ~\*]5 


2x 

a> X*/0 

R Reynolds number ( 

\ Mi® 


% Po*u 0 *8 

<%* 


M 0 

l Un ^ 

Mach number I , — 

\-yfvR*T 9 *, 

<x 

Prandtl number ==^ 


1. PRELIMINARY CONSIDERATIONS 


In the first phase of this investigation (reference 8) the 
stability of the laminar bound ary-layer flow of a gas is 
analyzed by the method of small perturbations, which was 
already so successfully utilized for the study of the stability 
of the laminar flow of an incompressible fluid. (See refer- 
ence 5.) By this method a nonsteady gas flow is investigated 
in which all physical quantities differ from their values in a 
given steady gas flow by small perturbations that are func- 
tions of the time and space coordinates. This nonsteady 
flow must satisfy the complete gas-dynamic equations of mo- 
tion and the same boundary conditions as the given steady 
flow. The question is whether the nonsteady flow damps to 
the steady flow, oscillates about it, or diverges from it with 
time — that is, whether the small perturbations are damped, 
neutral, or self-excited disturbances in time, and thus whether 
the given steady gas flow is stable or unstable. The analysis 
is particularly concerned with the conditions for the existence 
of neutral disturbances, which mark the transition from stable 
to unstable flow and define the minimum critical Reynolds 
number. 
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In order to bring out some of the principal features of the 
stability problem without becoming involved in hopeless 
mathematical complications, the solid boundary is taken as 
two-dimensional and of negligible c urvature and the boundary- 
layer flow is regarded as plane and essentially parallel ; 
that is, the velocity component in the direction normal to 
the surface is negligible and the velocity component parallel 
to the surface is a function mainly of the distance normal 
to the surface. The small disturbances, which are also two- 
dimensional, are analyzed into Fourier components, or normal 
modes, periodic in the direction of the free stream; and the 
amplitude of each one of these partial oscillations is a function 
of the distance normal to the solid surface, that is, 

u*'= Ua* 

In the study of the stability of the laminar boundary layer, 
it will be seen that only the local properties of the “parallel” 
flow are significant. To include the variation of the mean 
velocity in the direction of the free stream or the velocity 
component normal to the solid boundary in the problem 
would lead only to higher order terms in the differential 
equations governing the disturbances, since both of these 
factors are inversely proportional to the local Reynolds num- 
ber based on the boundary-layer thickness. (See, for ex- 
ample, reference 2.) By a careful analysis, Pretsch has 
shown that even with a pressure gradient in the direction of 
the free stream the local mean-velocity distribution alone 
determines the stability characteristics of the local boundary- 
layer flow at large Reynolds numbers (reference 9). Such 
a statement applies only to the stability of the flow unthin 
the boundary layer. For the interaction between the bound- 
ary layer and a main “external” supersonic flow, for example, 
it is obviously the variation in boundary-layer thickness and 
mean velocity along the surface that is significant. (See 
reference 10.) 

The aforementioned considerations also lead quite natu- 
rally to the study of individual partial oscillations of the 
form f(y)e taU ~ ct) , for which the differential equations of 
disturbance do not contain x and f explicitly. These partial 
oscillations are ideally suited for the study of instability, for 
in order to show that a flow is unstable it is unnecessary to 
consider the most general possible disturbance; in fact, the 
simplest will suffice. It is only necessary to show that a 
particular disturbance satisfying the equations of motion and 
the boundary conditions is seif-excited or, in this case, that 
the imaginary part of the complex phase velocity c is positive. 

In reference 8 the differential equations governing one 
normal mode of the disturbances in the laminar boundary 
layer of a gas were derived and studied verj thoroughly. 
The complete set of solutions of the disturbance equations 
was obtained and the physical boundary conditions that 
these solutions satisfy were investigated. It was found that 
the final relation between the values of c, a, and i? that 
determines the possible neutral disturbances (limits of 
stability) is of the same form in the compressible fluid as in 
the incompressible fluid, to a first approximation. The basis 
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for this result is the fact that for Reynolds numbers of the 
order of those encountered in most aerodynamic problems 
the temperature disturbances have only a negligible effect 
on those particular velocity solutions of the disturbance 
equations that depend primarily on the viscosity (viscous 
solutions). To a first approximation, these viscous solutions 
therefore do not depend directly on the heat conductivity 
and are of the same form as in the incompressible fluid, 
except that they involve the Reynolds number based on the 
kinematic viscosity near the solid boundary (where the vis- 
cous forces are important) rather than in the free stream. 
In this first approximation, the second viscosity coefficient, 
which is a measure of the dependence of the pressure on the 
rate of change of density, does not affect the stability of the 
laminar boundary layer. From these results it was inferred 
that at large Reynolds numbers the influence of the viscous 
forces on the stability is essentially the same as in an incom- 
pressible fluid. This inference is borne out by the results 
of the present paper. 

The influence of the inertial forces on the stability of the 
laminar boundary layer is reflected in the behavior of the 
asymptotic inviscid solutions of the disturbance equations, 
which are independent of Reynolds number in first approxi- 
mation. The results obtained in reference 8 show that the 
behavior of the inertial forces is dominated by the distribu- 
tion of the product of the mean density and mean vorticity 

p dy across fk e boundary layer. (The gradient of this quan- 
tity, or which plays the same role as the gradient 

of tRe vorticity in the case of an incompressible fluid, is a 
measure of the rate at which the ar-momentum of the thin 
layer of fluid near the critical layer (where w=c ) increases, 
or decreases, because of the transport of momentum by the 
disturbance.) In order to clarify the behavior of the inertial 
forces, the limiting case of an inviscid fluid (R— > => ) is studied 
in detail in reference 8. The following general criterions are 

obtained: (1) If the quantity vanishes forsome value 

of then neutral and self-excited subsonic disturb- 

ances exist and the inviscid compressible flow is unstable. 
(2) If the quantity does not vanish for some value 


then all subsonic disturbances of finite wave 


oiw > l ~W 

length are damped and the inviscid compressible flow is 
stable. (Outside the boundary layer, the relative velocity 
between the mean flow and the r-component of the phase 
velocity of a subsonic disturbance is less than the mean sonic 
velocity. The magnitude of such a disturbance dies out ex- 
ponentially with distance from the solid surface.) (3) In 
general, a disturbance gains energy from the mean flow 

if 


is positive at the critical layer (where w—c) and 
loses energy to the mean flow if { w=c 
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The general stability criterions for inviseid compressible 
flow give some insight into the effect of the inertial forces on 
the stability, but they cannot be taken over bodily to the 
real compressible fluid. Of course, if a flow is unstable in 
the limiting case of an infinite Reynolds number, the flow is 
unstable for a certain finite range of Reynolds number. A 
compressible flow that is stable when R— > <» , however, is not 
necessarily stable at all finite Reynolds numbers when the 
effect of viscosity is taken into account. One of the objects 
of the present paper is to settle this question. 

On the basis of the stability criterions ob tained in reference 8, 
some general statements were made concerning the effect of 
thermal conditions at the solid boundary on the stability of 
laminar boundary-layer flow. It is concluded from physical 
reasoning and a study of the equations of mean motion 


that the quantity 


vanishes for some value of 


if = 0; that is, if heat is added to the fluid through the 

solid surface or if the surface is insulated. If 


and is sufficiently large, that is, if heat is withdrawn from the 
fluid through the solid surface at a sufficient rate, the quan- 


tity ^y(yP^) never vanishes. Thus, when =0, the 

laminar boundary-layer flow is destabilized by the action 
of the inertial forces but stabilized through the increase of 


kinematic viscosity near the solid surface. When 


(D>°- 


the reverse is true. The question of which of these effects is 
predominant can be answered only by further study of the 
stability problem in a real compressible fluid. 

In the present paper this investigation is continued along 


the following lines: 


(1) A study is made of how the general criterions for 
instability in an inviseid compressible fluid are modified by 
the introduction of a small viscosity (stability at very large 
Reynolds numbers). 

(2) The conditions for the existence of neutral disturb- 
ances at large Reynolds numbers are examined (study of 
asymptotic form of relation between eigen- values of c, a, 
and R). 

(3) A relatively simple expression for the approximate 
value of the minimum critical Reynolds number is derived; 
this expression involves the local distribution of mean 
velocity and mean temperature across the boundary layer. 
This approximation will serve as a criterion from which the 
effect of the free-stream Mach number and thermal con- 
ditions at the solid surface on the stability of laminar 
boundary-layer flow is readily evaluated. The question of 


the relative influence of the kinematic viscosity and the 
distribution of p C j^ on stability would then be settled. 

(4) The energy balance for small disturbances in the real 
compressible fluid is considered in an attempt to clarify the 
physical basis for the instability of laminar gas flows. 

(5) In order to supplement the investigations outlined 
in the four preceding paragraphs, detailed calculations are 
made of the limits of stability, or the curve of a against 7? 
for the neutral disturbances for several representative cases 
of insulated and noninsulated surfaces. The results of the 
calculations are presented in figures 1 to 8 and tables I 
to IV. The method of computation of the stability limits is 
briefly outlined in reference 8, although the calculations 
were not carried out in that paper. 




Figure l.— Boundary-layer velocity profiles for insulated surface. Since <r is taken equal to 
unity, the temperature profile is given by jP=ITi— j^2 T i— Jig - - 2 * MqHp*, 
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Figure 2.- ^-Boundary-layer velocity profiles for nonlnsulated surface. M»=0.70. Tl is the 
ratio of surface temperature (deg abs.) to free-stream temperature (deg ahs.). Inflection 
more pronounced and farther out into fluid for Xi=1.25 than for insulated surface (Ti=1.10). 
No inflection for Ti=0.70, 0.80 , 0.90. 

In the present investigation the work of Heisenberg 
(reference 1) and Lin (reference 5) on the stability of a real 
incompressible fluid is naturally an indispensable guide. In 
fact, the methods utih'zed in the present study are analogous 
to those developed for an incompressible fluid. 

The present paper is concerned only with the subsonic 
disturbances. The amplitude of the subsonic disturbance 
dies out rapidly with distance from the solid boundary. In 
other words, the neutral subsonic disturbance is an “eigen- 
oscillation” confined mainly to the boundary layer and exists 
only for discrete eigen-values of c, a, and R that determine 
the limits of stability of laminar boundary-layer flow. Dis- 
turbances classified in reference 8 as neutral “supersonic,” 
that is, disturbances such that the relative velocity between 


the ar-component of the phase velocity of such a disturbance 
and the free-stream velocity is greater than the local mean 
sound speed in the free stream, are actually progressive sound 
waves that impinge obliquely on the boundary layer and 
are reflected with change of amplitude. For disturbances 
of this type the wave length and phase velocity are obviously 
completely arbitrary (eigen-values are continuous), and these 
disturbances have no significance for boundary-layer 
stability. 

When the free-stream velocity is supersonic (Mo>l), the 
subsonic boundary-layer disturbances must satisfy the re- 
quirement that v 0 * — c*<a 0 * or c>l — -jj (for M 0 <X e^O). 

Now, by analogy with the case of an incompressible fluid it 
is to be expected that for values of c greater than some 
critical value of c 0 , say, all subsonic disturbances are damped. 
Thus, when M 0 )>1, there is the possibility that, for certain 
mean velocity-temperature distributions across the boundary 
layer, neutral or self-excited disturbances satisfying the 
differential equations of motion, the boundary conditions, 

and also the physical requirement that c>l— jy- cannot be 

found. In that event, the laminar boundary flow is stable 
at all Reynolds numbers. This interesting possibility is 
investigated in the present paper. 

2. CALCULATION OF THE LIMITS OF STABILITY OF LAMINAR 
BOUNDARY LAYER IN A VISCOUS CONDUCTIVE GAS 

In order that the complete system of solutions of the differ- 
ential equations for the propagation of small disturbances 
in the laminar boundary layer shall satisfy the physical 
boundary conditions, the phase velocity must depend on the 
wave length, the Reynolds number, and the Mach number 
in a manner that is determined entirely by the local distribu- 
tion of mean velocity and mean temperature across the 
boundary layer. In other words, the only possible subsonic 
disturbances in the laminar boundary layer are those for 
which there exists a definite relation of the form (reference 8) 

c=c{a, R, M 0 2 ) (1) 

Since a, R, and Al a 2 are real quantities, the relation (1) is 
equivalent to the two relations 


C r =C T (a, R, Mo) 

(la) 

Ci~Ct(a, R, Mo 2 ) 

(lb) 


The curve a(a, R, M 0 2 ')= 0 (or a=a(R, Mo 2 )) for the neutral 
disturbances gives the limits of stability of the laminar 
boundary layer at a given value of the Mach number. From 
this curve can be determined the value of the Reynolds 
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Figuee 3.— Distribution of across boundary layer, as expresed by the function 

v(c)*=-^j F wllio11 appears directly in stability calculation. The 

value- of c=co at which (1— 2X) e=0.580 (X is small) is a measure of the stability of a given 
laminar boundary-layer flow. 

p, - ■ 6 Wa)l'' !i 

number below which disturbances of all wave lengths are 
damped and above which self-excited disturbances of certain 
wave lengths appear in a given laminar boundary-layer flow. 

In reference 8, it is shown that the relation (1) between the 
phase velocity and the wave length takes the following form: 

E(a, c, M 0 2 )=F(z) (2) 


In equation (2), F(z) is the Tietjens function (reference 11) 
defined by the relation 


m- 


where 


+ J 

' V 3/2 Hi/3 (1> 

03. . i 

§(*f) 3/2 ] 


z 

f~‘ 

l ») 3/2 

dr 


faRw/ \ 113 , v 

(M) 


(3) 


(4) 


and the quantity I? 1/3 (1> is the Hankel function of the first kind 
of order )■{. The prime denotes differentiation with respect 
to y. The function E(a, c, M 0 2 ), which depends only on the 



(b) Noninsulated surface. jVfo=0.70. 

Figure 3.— Concluded. 

asymptotic inviscid solutions <p t and (section 4 of ref erence’8) 

and not on the Reynolds number, is defined as follows: 


(: yi—y c )E(a,c,M 0 2 ) = 


<Pn <Pi 2 ' -¥P<Pl 2 
<P 2 \ ‘Pn' J r@<P 22 


where 


T \<pn' FlpW\' C(p n , 

Tx—MJc 1 vu FPvu 

Ti<p2\ Cipix , , „ 

l/ 0 2 c 2 ' P 1 ^ 22 


i3=ai/l— M 0 2 (l— C) 2 ' 

Vu=<Pi {Vi) 

i,j= 1,2 


(5) 


( 6 ) 


and yi and y 2 are the coordinates of the solid surface and the 
“edge” of the boundary layer, respectively. 

The Tietjens function was carefully recalculated in refer- 
ence 8, and the real and imaginary parts of the function 

$ (s) = y_ are plotted in figure 9. (The function $( 2 ) is 

found to be more suitable than F(z) for the actual calcula- 
tion of the stability limits.) 
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Figure i— Continued. 
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Re 


(fc) Insulated surface. Ma^O. Note that a-»as?^0 as i?-»oo for insulated surface. 
(1) Non insulated surface. jV/a=G.70. Note that a-tasj^Q as R-^<=a when Ti=1.25. 


Figuee 4. — Concluded. 
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The inviseid solutions tpi and cp 2 were obtained as power 
series in a 2 as follows (section 6 of reference 8) : 

p 1 (y;a i ,c,M 0 T ) = (w - c) (y;c, M 0 2 ) (7) 


<P2(y;a 2 ,c,M 0 2 ) = {w—c) ^a 2n k 2n +i (y;c,M 0 -) (8) 

71=0 


where for n > 1 




( W — C )' 


h 2n _i(y;c,M Q 2 )dy 


h a =1.0 


O 12 3-45 

M 0 

Figi. ee 8.— Critical temperature ratio Ti„ for stability of laminar boundary layer ajainst 

Mach number Aft. 

c.6 .e 


22 .< 


*r(*l *<(*} 


1.2 J 


1.0 c 


8 -./ 



20 2.4 2.8 3.2 3.6 4.0 4.4 4.8 


Fighbe 9— The functions $,( 2 ) and 4, ( 2 ). 


and for n> 1 


k 2u+1 (y;c,M a 2 )~ J] [ (w _ c)2 

Mo 2 ] dy J] ^^h^(y;c,Mo 2 )dy (10) 

and _ 

k 1 (y;c,M 0 2 ) = J] \_(w-cy~ M °*\ dy 

The lower limit in the integrals is taken at the surface 
merely for convenience. When y^>y e , the path of integra- 
tion must be taken below the point y=y e in tbe complex 
y-plane. The power series in a 2 is then uniformly conver- 
gent for any finite value of a. ' 

At the surface, the inviseid solutions are readily evaluated 

<pn — — C 'j 

< Pn'=w 1' 


^l' = --^(2W/o 2 c 2 )J 

At the “edge” of the boundary layer, the inviseid solutions 
are most conveniently expressed as follows: 

*\ 

<?12=(1— c) S a 2 "H 2n ic, M 0 2 ) 

71 = 0 

- ®22=(1 — c) S a in K la+ i(c, M 0 2 ) 

*u'=(l-c) S Mo 2 ) ( 5 

<Pn'=(l—c) [ 1 ~fff c 1 ) 7 C)2 ]g 0 ^K 2n (c,M 0 2 ) 
where 

B 2 ,, (c ,Mq 2 ) = A 2 U ( 2 / 2 ; 0 ,M 0 -) 

So=1.0 

K 2n+ l (c,M 0 2 ) = A*2 jj + 1 ( 2/2 ; C ,Mo 2 ) 

g 2 n-i(c,xMo 2 ) =[ j ^z^7 £ ^ ]~ 1 W (y» ;c,M a 2 ) - (13) 

(C,M 0 2 ) = [ 1 --W(1- C)2 ]~ 1 i 2 , +1 / (y 2 ;C,M 0 2 ) 


hn+i 


Kn — 1.0 
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With the aid of equations (11), the expression for 
E ( a , c, M 2 ) can be rewritten as follows: 


E(a,C,Mq 2 ) = 


($22 ~h AP 22 ) 


1+X(c) 


T 

Wl / (^22 , + /3'p22)+~r (‘Pli Pvn) 


where 




(14) 


(15) 


The relation (2) between the phase velocity and the wave 
length is brought into a form more suitable for the calculation 
of the stability limits by making use of the fact that for real 
values of c the imaginary part of E (a, c^M 0 2 ) is contributed 
largely by the integral K x {c, M/) . (The procedure to be 
followed is identical with that used by Lin in the limiting 
case of the incompressible fluid (reference 5, part III).) 
Define the function $( 2 ) by the relation 


$(z) = 


1 -m 


(16) 


Then, 

where 


$( 2 ) = 


(u+iv) 


WE-< 1+> > 1WR5 


(swfS) 


Equation (17) is equivalent to the two real relations 

(1+X)r ■ 


$»(s) = 


4> 


( 2 )— (l + X) ■ 


(l + Xu) 2 +XV 

u(l-f Xu) -{-Xr 2 " 
(l + Xu) 2 +XVj 


(17) 

(18) 

(19) 

( 20 ) 


The real and imaginary parts of $( 2 ) are plotted against 2 in 
figure 9. 

The dominant term in the imaginary part of the right-hand 
side of equation (18), which involves K x (c, M 0 2 ), is extracted 
by means of straightforward algebraic transformations. Re- 
lation (18) becomes 


u+iv= 


Wi c 

TV 


K, 




VI— M 0 2 (l—c) 2 
(1— c) 2 


/ CO \ CD 

1-Sa%-S 

\ n = 1 / n = 1 


a 2n+ W 2n+I 




(aH 1 -±a^A^ +1 ') 


( 21 ) 


where 

and for n=z 3 
and 

When c is real, 


N 2 =H 2 

N n =KiH n ^—K n 
M n =H 2 H n - 2 -H n " 

v~^l. V. Kl 


(22a) 

(22b) 


for those values of a and c that occur in the stability calculations. (This approximation is justified later in appendix A,) The 
imaginary part of the integral Ki{c, M a 2 ) is readily computed. It is found that 


Tc_ (wl_T/\ 

*(«>.')* W TJ 


(23) 


Now X(c) is generally quite small, therefore $j(z) can be taken equal to v(c ) and ^ r (z) can be taken equal to u as a 
zeroth approximation. From equations (19) and (20), when e.is real, 


1t l0) = $ r <#) (2 (0) ) 


(24) 

(25) 


By equation (24), z (0> is related to c with the aid of figure 9; and by equation (25), u (0) is also related to c. The quantity ceR 
is connected with c by means of the identity 


aR = 


Vc ( ZWi 

'w/(l + X) 3 \ c 


)‘ 


(26) 


and the corresponding values of a are obtained from equation (21) (slightly transformed) by a method of successive 
approximations. Thus 


( 27 ) 


f (1- 

~a 2 H,) 

-Vl-Af 0 *(l-c) a 

( 1 — c) 2 

fi-±J 

\ n=I / 

-s « 2n+l N 2a+1 

71 = 1 _ 

(u—L) 

•Wi 

1 

1 — 1 

1 1 

a. 2n M ! ^ 4f 0 2 (l c) 2 

“ ■ L ' 1 2n 1 c ^2 

^ alii - 

-S a 2 ^M. in+l ) 

71 = 1 / 
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where 



(The symbols M k and N t now designate the real parts of the 
integrals and A 7 *.) The iteration process is begun by 
taking a suitable initial value of a on the right-hand side of 
equation (27). The methods adopted for computing these 
integrals 'when the mean velocity-temperature profile is 
known are described in appendixes A to C. 

For greater accuracy, the values of 2 and u for a given real 
value of c are computed by successive approximations. From 
equations (19) and (20), 




(1 + X )p 

(H-\u w ) 2 -i-\V 


(28) 


,(»+!> 


= < j, r (,+t>( 2 C=-H)) 


(i+Xxt w ) 2 +xV 
_(1 + X)(l+Xtt‘*>)_' 


\v 2 


l+Xu 1 ’ 0 


(29) 


The value of v is always approximated by relation (24). 

Curves of wave number against Reynolds number for the 
neutral disturbance have been calculated for 10 representa- 
tive cases (fig. 4), that is, insulated surface at Mach numbers 
of 0, 0.50, 0.70, 0.90, 1.10, and 1.30 and heat transfer across 
the solid surface at a Mach number of 0.70 with values of the 
ratio of surface temperature to free-stream temperature T x of 
0.70, 0.80, 0.90, and 1.25. (It is found more desirable to base 
the nondimensional wave number and the Reynolds number 
on the momentum thickness 6, which is a direct measure of 
the skin friction, rather than on the boundary-layer thickness 
5, which is somewhat indefinite.) 

In figure 5 the minimum critical Reynolds number 
Rs , or the stability limit, is plotted against Mach num- 

ber for the insulated surface: and in figure 6(a) Re„ is 

plotted against Ti for the cooled or heated surface at a Mach 
number of 0.70. The marked stabilizing influence of a with- 
drawal of heat from the fluid is clearly evident. Discussion 
of the physical significance of these numerical results is 
reserved until after general criterions for the stability of the 
laminar boundary layer have been obtained. 


3. DESTABILIZING INFLUENCE OF VISCOSITY AT VERY LARGE 
REYNOLDS NUMBERS; EXTENSION OF HEISENBERG’S 
CRITERION TO THE COMPRESSIBLE FLUID 

The numerical calculation of the limits of stability for 
several particular cases gives some indication of the effects 
of free-stream Mach number and thermal conditions at the 
solid surface on the stability of the laminar boundary layer. 
It would be very desirable, however, to establish general 
criterions for laminar instability. For the incompressible 
fluid, Heisenberg has shown that the influence of viscosity is 
generally destabilizing at very large Reynolds numbers 
(reference 1). His criterion can be stated as follows: If a 

neutral disturbance of nonvanishing phase velocity and 
finite wave length exists in an inviseid fluid (R -><*>) for a 
given mean velocity distribution, a disturbance of the 
same ware length is unstable, or self-excited, in the real fluid 
at very large (but finite) Reynolds numbers. 

The same conclusion can be drawn from PrandtTs discussion 


of the energy balance for small disturbances in the laminar 
boundary layer (reference 12). 

Heisenberg’s criterion is established for subsonic disturb- 
ances in the laminar boundary layer of a compressible fluid 
by an argument quite similar to that which he gave originally 
for the incompressible fluid and which was later supple- 
mented by Lin (reference 5, part III). At very large 
Reynolds numbers, the relation (1) between the phase 
velocity and the wave length can be considerably simplified. 
When X is finite and c does not vanish, |^[ 1 at large 

Reynolds numbers. The asymptotic behavior of the Tietjens 
function F{z) as [zj— s-cn is given by (reference 5, part I) 


(Vi -Vc)F{z) = 


and the relation (1) becomes 


rw 

\ a ~ c 
V v c 


(yi—y e )E(a, c, M 0 2 )=E 1 (a, c, A/ 0 2 ) : 


— e 7174 


R 

a/ a ~ c 

V V; 


(30) 


(31) 


where E(a, c, fi/ 0 2 ) is given by equation (14). 

Suppose that a neutral disturbance of nonvanishing 

9 ? r 

wave number a s = : w an d phase velocity c t )>l — xr exists in 

A* ~ SXl c 

the inviseid fluid (limiting case of an infinite Reynolds 
number). The phase velocity c is a continuous function of R, 
and for a disturbance-of given wave number a, the value of c 
at very large Reynolds numbers will differ from c t by a small 
increment Ac. Both sides of equation (31) can be developed 
in a Taylor’s series in Ac, and an expression for Ac can be 
obtained as follows: 

E l (a,c,M<> 2 )=E l (cc s ,c 1 ,M a 2 )+ ^Ac+ . . . 



[1 + 0 (Ac)] 


The boundary condition 


(32) 


<p 22 r (a s , c„ 14 ! ) + ft? 22 ( ffi i; c„ il/ 0 2 ) = 0 (33) 

must be satisfied for the inviseid neutral disturbance, and 
the function E, (a„ c», M 0 2 ) vanishes (equation (14)) . Recog- 
nizing that 



reduces equation (32) for Ac to the form 


An 


From equation (14), 



(34) 



Cj 2 [<Pn (a», C, M a *)+P<pzi(a„ C, A/ 0 *>]J e _ e 

T\ ‘P12 (a*, c } , M/) + P s <p l2 (a s , C s , M 0 2 ) 


( 35 ) 
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By equations (12) and the boundary condition (33), the quantity (^r) c a is evaluated as follows: 


- n=l n=l yl — illp ( 1 — C s ) 

(1 -c s ) 2 fj a* n+t H in {c s ,Mo) + Vl-Wl-C.) 1 S l/o-) 

n=0 n=I 


fdE A _ 

\ 2>c J c s, a $ T\ 




( 36 ) 


where the primes now denote differentiation with respect to c. For small values of c s and a 3) the quantity ( ^f) c a ' s given 
approximately by the relation ... . . . .. . .. ... . 



cl 

Elo 


2-M 0 2 (i -c s y 


r (i— c s ) s vi— 4/ 0 2 (i — c s y 


-f-K/(c s 


Mo 2 ) 


(37) 


and the expression for Ac is 


Ac— 


T x 


a s -\ll — M 0 2 (l . — c,)V rt/4 


B ,2-iW(l-c s ) 2 

M as vJ s a-csf 


s Vl — -4Fg 2 (l —Cs) 2 Ki'(c„ Mq 2 ) 


Evaluation of the integral K x {c, A/ 0 2 ) yields the following result: 

7) , T t 


Ki (c, Mo 2 ) ■- 


u\'c (w, 


T Jy G 1 ). *=. (ln c ^ )+0(1) 


Since the quantity (vp 


vanishes (reference 8) , differentiation of equation (39) gives 


K ''< c -cW)=^+(|{^[|(T)]„,,j),.h lnCl -" >+0<1) 


(38) 


(39) 


(40-X 


Thus, Ki'(c„ Ml) is approximately real and positive for 
small values of c s . With c,>l— 4- ?• Ac must also be 

positive (equation (38)); therefore, a subsonic disturbance 
of wave length A^O, which is neutral in the inviscid com- 
pressible fluid, is self-excited in the real compressible fluid 
at very large (but finite) Reynolds numbers. 

In reference 8, it was proved that a neutral subsonic 
boundary-layer disturbance of nonvanishing phase velocity 
and finite wave length exists .in an inviscid compressible fluid 
d f dw\ 

\ P Ty)™ 


only if the quantity 


wT>l - 


M, 


dy 


vanishes for some value . of 


.-- If this condition is satisfied, then self-excited 


subsonic disturbances also exist in the fluid, and the laminar 
boundary layer is unstable in the limiting case of an infinite 
Reynolds number. By the extension of Heisenberg’s 
criterion to the compressible fluid, it can be seen that, far 
from stabilizing the flow, the small viscosity hi the real fluid 
has, on the contrary, a destabilizing influence at very large 
Reynolds numbers. Thus, any laminar boundary-layer 
flow in a viscous conductive gas for which the quantity 


(p My) van ^ ies f° r some value of w> 1 — ^ is unstable at 


d_ 
dy 

sufficiently high (but finite) Reynolds numbers, 
d f dw 


Unless the condition. ^ yp — 0 for some value w> I — 

is satisfied, all subsonic disturbances of finite wave length are 
damped in the limiting ease of infinite Reynolds number, and 
the inviscid flow is stable. Since the effect of viscosity is 
destabilizing at very large Reynolds numbers, however, a 
laminar boundary flow that is stabic.in the limit of infinite 
Reynoldanumber is not necessarily stable at large Reynolds 
numbers when the viscosity of the fluid is considered. 
(See fig. 4 (1).) In fact, for the incompressible fluid, Lin 
has shown that, every laminar boundary-layer flow is urn 
stable at. sufficiently high Reynolds numbers, whether or not 

tlxe vorticity gradient vanishes (reference 5, part III). 

In order to settle this question for the compressible fluid in 
general terms, the relation (1) between the complex phase 
velocity and the wave length at large Reynolds numbers 

must now be studied for flows in which the quantity 


i_ ( dvd\ 
dy V dy) 


does not v: 


vanish for any value of w>l — 

IVIq 
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4. STABILITY OF LAMINAR BOUNDARY LAYER AT LARGE 
REYNOLDS NUMBERS 

The neutral subsonic disturbance marks a possible 
“boundary” between the damped and the self-excited 
disturbance, that is, between stable and unstable flow. 
Thus, the general conditions under which self-excited dis- 
turbances exist in the laminar boundary layer at large 
Reynolds numbers can be determined from a study of the 
behavior of the curve of a against R for the neutral disturb- 
ances. When the mean free-stream velocity is subsonic 
Oio<l) , the physical situation for the subsonic disturbances 
at large Reynolds numbers is quite similar to the analogous 
situation for the incompressible fluid. The curve of a 
against R for the neutral disturbances can be expected 
to have two distinct asymptotic branches that enclose a 
region of instability in the aR-plane, regardless of the local 
distribution of mean velocity and mean temperature across 
the boundary layer. When the mean free-stream velocity 
is supersonic (V„>1), the situation is somewhat different; 
under certain conditions (soon to be defined) a neutral or 

a self-excited subsonic disturbance ( 1 A) cannot exist 

at any value of the Reynolds number. For this reason, it 
is more convenient to study the case of subsonic and super- 
sonic free-stream velocity separately. 


A. SUBSONIC FREE-STREAM VELOCITY (Me<D 


The asymptotic behavior at large Reynolds numbers of 
the curve of a against R for the neutral disturbances is 
determined by the relations (19) to (22) between a, R, and c 
for real values of c. For small values of a and c, these rela- 
tions are given approximately by 


i’(c)=$i(z) 


ttWiC T c 2 r d 
T x (w/Y \_dy V TJX= C 


(41) 


U = ^ r (z) 


(42) 


r T 1.76 / -A 3 

E= a( I c) {Wl,y (43) 


a= TT (44) 

As R — * a> , either z— > or 2 remains finite while both a and 
c approach 0. These two possibilities correspond to two 
asymptotic branches of the curve of a against R. 

lower branch. — If z remains finite as !?—*<», then c-^0; 
and by equation (41), $j(z)— M3. Therefore, 2 — *2.2-9 while 
u— *2.29 (fig. 9). From equations (43) and (44), along the 
lower branch of the curve of a against R for neutral stability 


Upper branch. — Along the upper branch of the curve of a 
against R for neutral stability, and 


*i(z) 


XlL\'c 

T/ 

~d 

1 

~ 2i 

(w/Y 

_dy V TJ_ 

w=c y2? 


while u-*1.0 (fig. 9 and equation (42)). If the quantity 

54 (} ~T ) c ^ oes n °t Yan i s i L f° r any value of tt’)> 0, then by 

equation (47) c must approach zero as 2 — *co. Along this 
branch, 


R ■■ 


«) u (1- 

-Wr 

-2 ]_ 

2rTy- 2i )T d 

{ w '\~ 
\TJ_ 

| 2 cr 6 

t) 

^ 



tt’Vv'i— - 

w a 



(48) 


(49) 


and a— *0 at large Reynolds numbers (fig. 4 (I)). 


On the other hand, if ^ \Tj can ^ es f° r some value of 

u’=c,)> 0, then by equation (47) c— >c s and a— *a, as both z and 
R approach . Now, 

r d /w'm jr&_ vzo * i _ 

Idy l. T)\ a= rw V t) 1 w/ T \ [dtf \ T)i 

rd 2 W:" ) c’—c> , ,, n , 

to/ ( 2(iV)* 4 '*’ (50) 

“&(?)} does not vanish, (see appendix D), then by 

equations (44) and (47), along the upper branch of the curve 
of a against R for the neutral disturbances, 


R-- 


(w x 'Y 


1 _ 


'2 v 2 !?- 24 ( fii (yy\ 

Iwwi 


\ s aC s (c-C t ) 2 

) 


a^^r- y'l— V 0 2 (l— c) 2 


(51) 


(52) 


and c— *c s ^0, a-*a 3 ^0 at large Reynolds numbers (figs. 4 (k) 


and 4 (1)). If 
replaced by 

2(w l ') 10 


Vdd_ /v/\~ 

ldy 2 \T) Ji 


vanishes, the relation (51) is 


R-- 


1 


1 


VT?- 2 ' Vd) (w'X 
\W\T)_ 


) 2 aC 5 (c 2 — C s 2 ) 2 


(53) 


R-- 


(w/ra-AiY ) 3 ' 2 1 


TV 


c — 2.29 — ~ 


T, 


Wl-A/ 0 2 

and a— >0 at large Reynolds numbers (fig. 4 (1)). 


(45) 

(46) 


which reduces to the relation obtained by Lin in the limiting 
case of an incompressible fluid when A/ 0 -*0, the solid bound- 
ary is insulated, and w"~Q for some value of w=c { )>0. 
(See equation (12.22) of reference 5, part III.) 


If the quantity 


d_ 

dy 



vanishes at the solid boundary 
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(that is, for w=0) , it can be shown from the equations of 


motion (appendix D) that 


d 2 / yfff 


is always positive— 


except in the limiting case of an incompressible fluid. For 
small values of y, the quantities 


and are both 


For large values of y, however, 
w' 


positive and increasing. 

-yr-i»0, physically; therefore ^ must have a maximum, or 

d 

dy V T 


= 0 for some value of 0, and this case is no 

different from the general case treated in the preceding 
paragraph. In the limiting case of an incompressible fluid, 

when w" vanishes at the surface, since w/" 

always vanishes in this case. From equation (48) the rela- 
tion between a and R along the upper branch of the neutral 
stability curve is therefore 


R-- 


2 (O 1 


1 


1 


(54) 


7T 2 (W] lv ) 2 a 10 

which is identical with equation (12.19) in reference 5, part III. 

Thus, regardless of the behavior of the quantity ^ 

regardless of the local distribution of mean velocity and mean 
temperature across the boundary layer — when M 0 <1, the 
curve of a against R for the neutral disturbances has two 
distinct branches at large Reynolds numbers. From physi- 
cal considerations, all subsonic disturbances must be damped 
when the wave length is sufficiently small (a large) or the 
Reynolds number is sufficiently low. Consequently, the 
two branches of the curve of a against R for the neutral dis- 
turbances must join eventually, and the region between 
them in the off -plane is a region of instability; that is, at a 
given value of the Reynolds number, subsonic disturbances 
with wave lengths lying between two critical values \i and 
X 2 (m and a 2 ) are self-excited. Thus, when M 0 <{1, any lami- 
nar boundary-layer jlow in a viscous conductive gas is unstable 
at sufficiently high {but finite) Reynolds numbers. 

The lower branch of the curve, of a against R for the 
neutral disturbances is virtually unaffected by the distribu- 

. f d ( w ' 
iron of 


across the boundary layer, but for the upper 


branch the behavior of the quantity 


d_(vfi 

dy\T 




is decisive. 


When 


ly ^ or some va ^ ue of w=c 5 >0, the neutral 

subsonic disturbance passes continuously into the charac- 
teristic inviscid disturbance c=c s and a=a s as. I?— > co . 
This result is in accordance with the results obtained in 
reference 9 for the inviscid compressible fluid and is in agree- 
ment with Heisenberg’s criterion. In addition, all subsonic 

2ir 

disturbances of finite wave length Xj>X s =— (and nonvanish- 


ly Ct) d° es n °t van ‘ s ^ f° r an y valoe of 0, then except for 

the “singular” neutral disturbance of zero phase velocity 
and infinite wave length (c=0 and «=0), all disturbances 
are damped in the inviscid compressible fluid. This 
singular neutral disturbance can be regarded as the limiting 
case of the neutral subsonic disturbance in a real compressi- 
ble fluid as 2?— » “ . 

B. SUPERSONIC FREE-STREAM VELOCITY (M»>I) 

When the velocity of the free stream is supersonic, the 
subsonic boundary-layer disturbances must satisfy not only 
the differential equations and the boundary conditions of the 

problem but also the physical requirement that e r )> 1 — rr 

nr o 

The asymptotic behavior at large Reynolds numbers of the 
curve of a against R for the neutral subsonic disturbances 
is determined by the approximate relations (41) to (44), with 

the additional restriction that c)>l— w- As c— »1 — 

AIq iuo 

0 by equation (44); therefore, !?->«> by equation (43). 

The corresponding value (or values) of z is determined by 

equation (41) as follows: 


$ t .( 2 )=y(c)=i;^l-ij 


r T 2 d (w'W 

““ (55) 

Now from physical considerations, (fifj < 0 for large 
values of y. Therefore, if ^ (yr)~ 0 (changes sign) for some 
value of w=c s >l— then, in general, j ~_dy(fip) \ 

10-1 Ms 

and $t(z) i <0 (equation (55)). From figure 9 it can bo 

C =1 7-y - 

Mo 

seen that in this case there is only one value of z {z u say) 
corresponding to the value of 4p(2) given by equation (55). 
From equations (42) to (44) along the lower branch of the 
curve of a against R for the neutral disturbances, 


R-. 


Wi 1 — 


'Mr 


27- 76 (w/)V 1 

0 -iJ ; 

V®, 


TxUi 


yrv 


l 

'Mg 


(56) 


(57) 


and c — >1 ■ 


'Mg 


at large Reynolds numbers (fig. 4(k) ) . The 


upper branch of the curve in tills case is given by equations 


ing phase velocity 0<c,<c s ) are self-excited in the limiting , 
case of infinite Reynolds number. On the other hand, when j A ams es ' W1 * ; c 


(51) and (52), or by equations (53) and (52.) if 

1 


~cP fw'\~ 


=>c s > 1 — and «- 

JSJ.Q 


>a,p^0. 
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If 


dy 


/w / \ 

\TJ 


vanishes for w—1 — yt' then 
Ai 0 


as R- 


along the upper branch of the curve of a against R for the 


neutral disturbances, and 3>j(z)- 


W{ 


<2 a - c s 
Vc 


Now a — >0 as 


1 


c — >1 — tt in this case also (equation (57) with ■w 1 =1.0) so 
AI 0 

that 




2(w-.') u M 0 - 

_2T 4.24 
V 1 1 


1 


1 - 


M n 


&(?)],!' 


(58) 


Along the lower branch of the cuiwe of a against R at large 
Reynolds numbers, a, R, and c are connected by equations 
(56) and (57), with z x = 2.29 and-u : =2.29. In spite of the 


fact that 


1— 171 a neutral sonic disturb- 
er 0 


of finite wave length does not exist- in 


dy (tf ) = 0 for w 
ance (c= 1 -^) 

the inviscid fluid unless R)(c) = |* 2 ~ Af & 2 J dy is 

positive. (See section 10 of reference 8 .) Calculation 
shows that K x (c) is almost always negative (equation (40)); 
therefore, in general, the sonic disturbance of infinite wave 
length (a= 0 ) with constant phase across the boundary 
layer exists only in the inviscid fluid (R — *■<=). 

If i^) does n °k van i s ^ f° r an 7 value of wA 1 — it 
is certain that 1 <0 and by equation (55) 




®‘C g )«- 1 -A>° 

Ale 


When v x_<0.580 (approx.), there are two 

-Us 

values of 2 (z a and Zj, say, with 2 3 >z 2 ) corresponding to the 
value of $i(s) given by equation (55). (See fig. 9.) Along 
the two asymptotic branches of the curve of a against R for 
the neutral disturbances, a, R, and c are connected by rela- 
tions of the form of equations (56) and (57), with 2 and u 
replaced by z 2 and u 2) respectively, along the lower branch 
and by z 3 and u 3 , respectively, along the upper branch. At 
a given value of the Mach number, the value of v i_ is con- 

1 Mi 

trolled by the thermal conditions at the solid surface. (See 

seetion 6 .) When these conditions are such that v i_«0.580, 

Ml 

then £2 = 23 , and the two asymptotic branches of the 
curve of a against R for the neutral disturbances coincide. 

When v j_A0.580 (approx.), it is impossible for a neutral 
1 Ml 

or a self-excited subsonic disturbance to exist in the laminar 
boundary layer of a viscous conductive gas at any value of 


the Revnolds ntimber. In other words, if 


L 


AO. 580 


(approx.), the laminar boundary layer is stable at all values 
of the Reynolds number. (Of course, in any given case, the 
critical conditions beyond which only damped subsonic dis- 
turbances exist can be calculated more accurately from the 
relations (28) and (29). See section 5 on minimum critical 
Reynolds number.) 

The preceding conclusion can also be deduced, at least 
qualitatively, from the results of a study of the energy 


balance for a neutral subsonic disturbance in the laminar 
boundary layer. A neutral subsonic disturbance can exist 
only when the destabilizing effect of 'viscosity near the solid 
surface, the damping effect of viscosity in the fluid, and the 
energy transfer between mean flow and disturbance in the 
vicinity of the inner “critical layer” all balance out to give 
a zero (average) net rate of change of the energy of the dis- 
turbance. (See Schliehting’s discussion for incompressible 
fluid in reference 4.) In reference 8 it is shown that the sign 
and magnitude of the phase shift in u* r through the inner 
“critical layer” at w=c is determined by the sign and magni- 
tude of the quantity • The corresponding 


apparent shear stress r c *= —p* which is zero for w<Cc 

in the inviscid compressible fluid, is given by the following 
expression for wT>c (reference 8 ): 


ur=Pa(« 


2 a 
2 


( 5) 3 [dy ( 


(59) 


If the quantity |j^ (vp)], r=i is negative, the mean flow ab- 


sorbs energy from the disturbance 




is positive, 


energy passes from the mean flow to the disturbance. In 
the real compressible fluid, the thickness of the inner critical 
layer in which the viscous forces are important is of the order 

j^ T/ 5 ) an d the phase shift in u*' is actually brought about 


(■ 


“rj 


by the effects of viscous diffusion (^of the quantity p 
through this layer. ■ 

As shown by Prandtl (reference 12 ), the destabilizing ef- 
fect of viscosity near the solid surface is to shift the phase of 
the “frictional” component u fT *' of the disturbance velocity 
against the phase of the “frictionless” or “inviscid” compo- 
nent u in ,*' in a thin layer of fluid of thickness of the order of 

■ — By continuity, the associated normal component 


/-T 5 . (It was 


Vi 

v fr *' is of the order of hjpd / — ^ 1 

V ac— \j c- 

I Vl V Vi 

shown in section 1 of reference S that for large values of <xR 

the “frictional” components of the disturbance also satisfy 

the continuity relation ^ + g— ? = 0 in the compress- 
ible fluid.) The corresponding apparent shear stress 
n*= — pi*u*'v*' is given by the expression 


Ti —Po’ 


p 0* LV ■ uo * /ij ■ 


V 


CL 

~R. 

Ti 


(60) 


But from equations (11) 


IV uo* J 1 


\ 


Ti 

X 

! 

2 \— l/ 0 2 c 2 ® 21 


Tx 

c 


( 61 ) 
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n* ~ Pt>*{v>a *) 2 J j ft 
Mo- 


rn 


Since the shear stress associated with the destabilizing 
effect of viscosity near the solid surface and the shear stress 
near the critical layer act roughly throughout the same 
region of the fluid, the ratio of the rates of energy trans- 

(approximately |* t by the two physical 


ferred 
processes, is 


AM 


r * 

E l *\ 


Tl* 


7 TWiC 1 c 

”2 Ti ( w c J 


IF d /w' 


U3/2 


where 


\v(c)\z ™ 


(63) 


Z a ~ 


R 


V, M 2 

If the quantity ^ (^t) ne 9 a ^ ve an d sufficiently large 


when 


w— Ci, sajq then the rate at which energy is absorbed by the 
mean flow near the inner “critical layer” plus the rate at 
which the energy of the disturbance is dissipated by viscous 
action more than counterbalances the rate at which energy 
passes from the mean flow to the disturbance because of the 
destabilizing effect of viscosity near the solid surface. Con- 
sequently, a neutral subsonic disturbance with the phase 
velocity does not exist; in fact, all subsonic disturb- 
ances for which cScj are damped. 'When JV/ 0 <1, there is 
always a range of values of phase velocity O^c^c 0 for 
\E * 

which the ratio , given by equation (63), is small enough 

for neutral (and self-excited) subsonic disturbances to exist 
for Reynolds numbers greater than a certain critical value. 
When ,l/ : >] , however, because of the physical requirement 

thatc)>l— the possibility exists that for certain ther- 

0 _ 

mal conditions at the solid surface the quantity K (?)],.. 

\E *1 

is always sufficiently large negatively (and therefore is 

\E I 

sufficiently large) so that only damped subsonic disturbances 
exist at all Reynolds numbers. Of course, if (^t) van ishes 

for some value of w^l- it is certain that a(c)<(0.580 

IVIq 

for some range of values of the phase velocity 1— -igcg'co- 

1\1q 

In that case, neutral and self-excited subsonic disturbances 
always exist for R^>R CTm . and the flow is always unstable at 
sufficiently high Reynolds numbers, in accordance with 
Heisenberg’s criterion as extended to the compressible fluid 
(section 2). 


A discussion of the significance of these results is reserved 
for a later section, (section 6) in which the behavior of the 

quantity ^ (^t') be r(l lated directly to the thermal 

conditions at the solid surface and the free-slream Mach 
number. 

5. CRITERION FOR THE MINIMUM CRITICAL REYNOLDS 

NUMBER 

The object of the stability analysis is not only to deter- 
mine the general conditions under which the laminar bound- 
ary layer is unstable at sufficiently high Reynolds numbers 
but also i if possible, to obtain some simple criterion for the 
limit of stability of the flow (minimum critical Reynolds 
number) in terms of tbe local distribution of mean velocity 
and mean temperature across the boundary layer. For plane 
Couette motion (linear velocity profile) and plane Poiseuillc 
motion (parabolic velocity profile) in an incompressible fluid, 
Symge (reference 13) was able to prove rigorously that a 
minimum critical Reynolds number actually exists below 
which the flow is stable. His proof applies also to the 
laminar boundary 7 layer in an incompressible fluid, with only 
a slight modification (reference 5, part HI). Such a proof is 
more difficult to give for the laminar boundary layer in a 
viscous conductive gas; however, the existence, in general, 
of a minimum critical Reynolds number can be inferred 
from purely physical considerations. A study of the energy 
balance for small disturbances in tbe laminar boundary layer 
shows that the ratio of the rate of viscous dissipation to the 
rate of energy transfer near the critical layer is 1/7? for a 
disturbance of given wave length while the energy transfer 
associated with the destabilizing action of viscosity near the 
solid surface bears the ratio 1/V7? to the energy transfer near 
the critical layer. Thus, the effects of viscous dissipation 
will predominate at sufficiently low Reynolds numbers and 
all subsonic disturbances will be clamped. The two distinct 
asymptotic branches of the curve of a against 7? for the 
neutral disturbances at large Reynolds numbers must join 
eventually (section 4) and the flow is stable for all Reynolds 
numbers less than a certain critical value. 

An estimate of the value of R cr . , which will serve as a 
stability criterion, is obtained by taking the phase velocity e 
to have the maximum possible value c 0 for a neutral subsonic 
disturbance; that is, for cj>c Q all subsonic disturbances are 
damped. This condition is very nearly equivalent to the con- 
dition that all be a minimum, which was employed by Lin for 
tbe ease of the incompressible fluid (p. 285 of reference 5, 
part III). The condition c=c 0 occurs when. <1 q(z) is a maxi- 
mum; that is, when 4>j(z) =0.580, z 0 =3.22, and $ f (2o)=l-4S 
(fig. 9). The corresponding value of c=c 0 can be calculated 
from the relations (19) to (22). Neglecting terms in X s 
(X is usually very small) and taking 14=1.50 gives 

$i(z) ~[1— 2X(c)]y(c) 


(04) 
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where 


and 


v(e) = — w 


Wi'c 

~Tf 


T 2 d_ 
(w'Y dy 



\(c) 


wi'(yc—vi) 

c 


1 


(65) 

( 66 ) 


It is only necessary to plot the quantity (1— 2\)» against c 
for a given laminar boundary-layer flow and find the value 
of c=c 0 for which (1 — 2X)s=0.5SG. The corresponding value 
of aR is determined from the relation 

ai?=[T(c 0 )r-- 6 (uV) 2 (|y (67) 

and this value of aR is very close to the m inim um value of 
aR. A rough estimate of the value of a for c=c 0 is given by 
the following relation (equation (27)): 


a «w I 'c 0 Vl--W(l-c 0 ) 2 (68) 


This estimated value of a is, in general, too small. The fol- 
lowing estimate of R C r mbl is obtained by making an approxi- 
mate allowance for this discrepancy and by taking round 
numbers: 


Rc 


25 [T(cq)F 6 < 

CoST-iI/ 0 2 (l-Co) 2 


(69) 


or 


250- 


S,.. »- 


I Uo* 

il: rSZ 


"V VC)*X* 


[^(Co)] 1 


76 

\<>v /i 


c o 4 Vl — (1 c 0 ) 3 


(70) 


For zero pressure gradient, the slope of the velocity profile 
at the surface is given very closely by (appendix B) 


and 


Therefore 


(*). 


/bw\ 
' Ti 

0.332 
= Ti 


iJM- 

V 


Re 


= 0.6667 


[T(c 0 )Y^ 


Ti c 0 \>l-M 0 \l-c 0 y 


(71) 


The expression (71) is useful as a rough criterion for the 
dependence of Re on the local distribution of mean 

velocity and mean temperature across the boundary layer. 

It is immediatelv evident that Re — *•<» when e tt — M— tt - 

" cr min IVaq 

When [(1— 2X) i>] _ i Si 0.580, the laminar boundary layer is 
stable at all values of the Reynolds number. (This condition 
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is an improvement on. the stability condition v _ j Si 0.580 

(approx.) stated in section 4.) 

In the following tables and in figures 5 and 6 (a) the esti- 
mated values of Re c , m!a given by equation (71) can be com- 
pared with the values of Ri crmia taken from the calculated 
curves of as against Re for the neutral disturbances. For 
the insulated surface, the values are 


Mt 

Cl 

rcc) 

R i CImm 
test.) 

Cfig.p 

0 

0.4130 

1.0000 

195 

150 

.50 

. 4400 

1.Q40S 

170 

136 

.70 

.4600 

1.07S2 

150 

126 

.90 

.4350 

1.1254 

129 

115 

1.10 

.5139 

1.1803 

109 

104 

1.30 

.5450 

1.2406 

92 

92 


For the noninsulated surface when Af o =0.70, the values are 


T, 

C9 

net) 

(est.) 

■Ri 

(fig. i) 

0.70 

0. 1872 

0.7712 

5377 

5150 

.SO 

.2619 

.3716 

1463 

1440 

.90 

.3394 

.9562 

524 

523 

1.25 

-5104 

L1449 

S9 

63 


The expression (71) for Ri„ mfri gives the correct order of 
magnitude and the proper variation of the stability limit 
with Mach number and with surface temperature at a given 
Mack number. 

The form of the criterion for the minimum critical Rey- 
nolds number (equation (71)) and the results of the detailed 
stability calculations for several representative cases (figs. 3 
and 4) show that the distribution of the product of the 

density and the vorticity P across the boundary layer 

largely deter min es the limits of stability of laminar boundary- 
layer flow. The fact that the “proper” Reynolds number 
that appears in the boundary-layer stability calculations is 
based on the kinematic viscosity at the inner critical layer 
(where the viscous forces are important) rather than in the 
free stream also enters the problem, but it amounts only to 
a numerical and not a qualitative change, when the usual 
Reynolds number based on free-stream kinematic viscosity 
is finally computed. Whether the value of Rs CTmt!i for a 
given laminar boundary-layer flow is larger or smaller than 
the value of Rt„ min for the Blasius flow, for example, is 

determined entirely by the distribution of p ~ across the 
boundary layer. If the quantity 4 ( p 44 is negative and 


dy V dy 

large near the solid surface so that the quantity (1— 2X)t?(c) 
reaches the value 0.580 when the value of c=c 0 is less than 
0.4186, the flow is relatively more stable than the Blasius 


883026—5' 


-21 
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of across the boundary layer, which remained open in 


flow. If the quantity ^ (p is positive near the solid 

surface, so that (1 — 2X)»(c) =0.580 when w (or e)>0.4186, 
the flow is relatively less stable than the Blasius flow. Thus, 
the question of the relative influence on Rt„ min of the kine- 
matic viscosity at the inner critical layer and the distribution 
dw 
d y 

the concluding discussions of reference 8, is now settled. 

The physical basis for the predominant influence on 

Re C r min of the distribution of p ^ across the boundary layer 

is to be found in a study of the energy balance for a subsonic 
boundary-layer disturbance (section 4). The distribution of 
dw 
dy 

velocity c 0 or the maximum possible distance of the inner 
critical layer from the solid surface for a neutral subsonic 
disturbance. The greater the distance of the inner critical 
layer from the solid surface, the greater (relatively) the rate 
of energy absorbed by the mean flow, from the disturbance 
in the vicinity of the critical layer (equations (61) and (62)). 
When c 0 is large, therefore, the. energy balance for a neutral 
subsonic disturbance is achieved only when the destabilizing 
action of viscosity near the solid surface is relatively large or, 

1 


p t— determines the maximum possible value of the phase 


in other words, when 


<*o ' 




:«c 0 3/2 is large and the Reynolds 


V 

\ -c g 

number R 0 , which is very nearly equal to R C r min i is correspond- 
ingly small. On the other hand, when c 0 is. small and the 
inner critical layer is close to the solid surface, the rate at 
which energy is absorbed from the disturbance near the 
critical layer is relatively small and the rate at which energy 
passes to the disturbance near the solid surface, which is of 

the order of — = ; is also relatively small for energy bal- 


ance; consequently R c . 


is large. 


6. PHYSICAL SIGNIFICANCE OF RESULTS OF STABILITY 

ANALYSIS 

A. GENERAL 

From the results obtained in the present paper and in 
reference 8, it is clear that the stability of the laminar bound- 
ary layer in a compressible fluid is governed by the action of 
both viscous and inertia forces. Just as in the. case of an 
incompressible fluid, the stability problem cannot be under- 
stood unless the viscosity of the fluid is taken into account. 
Thus, whether or not a laminar boundary-layer flow is un- 
stable in the inviseid compressible fluid (R— »«>), that is, 
whether or not the product of the density and the vorticity 

p ( i^ has an extremum for some value of «C> l—Jh-, there is 
ciy .Mo 

always some value of the Reynolds number R, r below 

*■ *- ^min 

which the effect of viscous dissipation predominates and the 
flow is stable. On the other hand, at very large Reynolds 
numbers the influence of viscosity is deslablizing. If the 


free-stream velocity is subsonic, any laminar boundary-layer 
flow is unstable at sufficiently high (but finite) Reynolds 
numbers, whether or not the flow is stable, in the inviseid 
fluid when only the inertia forces are considered,. 

The action of the inertia forces is more decisive for the 
stability of the laminar boundary layer if the free-stream 
velocity is supersonic. Because of the physical requirement 
that the relative phase velocity (1— c) of the boundary-layer 

disturbances must be subsonic, it follows that c>l — -ry)>0 

if i 0 


and the quantity 


Vd ( 

jiy v 


dw y 
dyJ_\ K= 


can be large enough nega- 


tively under certain conditions so that the stabilizing action 
of the inertia forces near the inner critical layer (where 
w=c^> 0) is not overcome by the destabilizing action of 
viscosity near the solid surface. In that case, undamped 
disturbances cannot exist in the fluid, and the flow is stable 
at all values of the Reynolds number. 

Regardless of the free-stream velocity, the distribution of 

the product of the density and the vorticity p-^ across the 

boundary layer determines the actual limit of stability, or the 
minimum critical Reynolds number, for laminar boundary- 
layer flow in a viscous conductive gas (equation (71)). 

Since the distribution of p across the boundary layer in 

turn is determined by the free-stream Mach number and 
the thermal conditions at the solid surface, the effect of these 
physical parameters on the stability of laminar boundary- 
layer flow is readily evaluated. 

B. EFFECT OF FREE-STREAM MACH NUMBER AND THERMAL CONDITIONS 
AT SOLID SURFACE ON STABILITY OF LAMINAR BOUNDARY LAYER 

The distribution of mean velocity and mean temperature 
(^and therefore of across the laminar boundary layer in 

a viscous conductive gas is strongly influenced by the fact 
that the viscosity of a gas increases with the temperature. 
(For most gases, p<xT m (m= 0.76 for air) over a fairly wide 
temperature range.) When heat is transferred to the fluid 
through the solid surface, the temperature and viscosity 
near the surface both decrease along the outward normal, and 
the fluid near the surface is more retarded by the viscous 
shear than the fluid farther out from the surface — as com- 
pared with the isothermal Blasius flow. The velocity profile 
therefore always possesses a point of inflection (where w" — 0) 
when he.at is added to the fluid through the solid surface, 
provided there is no pressure gradient in the direction of the 

n a- d / dw\ w" w'T' ,, ... 

mam flow, bmee d/y\liy)~T~~T~ > “ ie c l uai “ffy 


( \ p< ~3y van ^ lcs an< ^ P^j ^ ias an extremum at some point 


d_ 
dy 

in the fluid. On the other hand, if heat is withdrawn from 
the fluid through the solid surface, and ^ arc both posi- 
tive near the surface and the fluid near the surface is less 
retarded than the fluid farther out — as compared with the 
Blasius flow. The velocity profile is therefore more convex 
near the surface than the Blasius profile. 
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As pointed out in section 11 of reference 8, the influence, of 
the variable viscosity on the behavior of the product of the 

density and the vorticity p ^ can be seen directly from the 

equations of motion for the mean flow. When there is no 
pressure gradient in the direction of the main flow, the fluid 
acceleration vanishes at the solid surface, or 


and 


( d 


St* \ _r a f— 

dy*) l Ldy* V 1 


j du 

by*JJ 


=0 


(72) 


/ b-u* \ 1_ f bu* \ _ __m /d T*\ /du*\ , 

\by* 2 J~ J^\by*J\by*J~~T\^)Xby*J 1 ^ 

Thus, when heat is added to the fluid through the solid sur- 
face (Ti'<^ 0), J is positive, and the velocity profile is 

concave near the surface and possesses a point of inflection 
for some value of u>^> 0; when heat is withdrawn from the 

fluid (T/>0), \^* 2 j is negative, and the velocity profile 

is more convex near the surface than the Blasius profile. 

The behavior of the quantity 


d% * 

by* 


y*J 

From equation (73), in nondi- 


m-fl 


Ti 2 


Ti'wi' 


(74) 


is parallel to that of 
mensional form, 

Vd/ du?\~j fw'X 

L dy V ■ p dy ) 'Jr ~ [dy \TJ_ 

Differentiating the dynamic equations once and making use 
of the energy equation gives the following expression for 

[$(t) 1 (a PP endLx rj i : 


-1)(y-1)M 0 Ti 


2 OhX + 2 (m+ 1) W 


Thus, for zero pressure gradient , 

Now, if the surface is insulated, the quantity 


w 
2V Y 


Ti® 
(75) 

is always positive. 


vanishes, but 


d 1 fw'X] \_ n A d /w'\ 
_dy*\T )\ >0 and dy \ t) ■ 


and ^ both in- 


crease with distance from the solid surface. Since 


w 


>0 


far from the solid surface, ^ has a maximum and 
vanishes for some value of uq>0. If heat is added to the 
fluid through the solid surface (2Y<0), is already 


>0, the quan- 


positive at the surface; and since ^^2 

tity ^ vanishes at a point in the fluid which is farther 

from the surface than for an insulated boundary at the same 
Mach number (figs. 3 (a) and 3 (b)). Consequently, the value 
of C — Cf, for which the function 


(1 2\)v(c) *(1 2X) ^ [(«,')* ( r)_ 

reaches the value 0.580 is larger than the value for the insulated 
surface. By equation (71), the effect of adding heat to the 
fluid through the solid surface is to reduce Re and to 

CT min 

destabilize the flow, as compared with the flow over an in- 
sulated surface at the same Mach number (fig. 6) . 

If heat is withdrawn from the fluid through the solid sur- 
face, Ti'YO and X ~( w 


\dy\T 
of heat transfer is 


01 is 


is negative. In fact, if the rate 
sufficiently large, the quantity 

dy (t) does n °t Yanish within the boundary layer (fig. 3 (b)). 

The value of c=c 0 for which the function (1— 2X)r(e) 
reaches the value 0.580 is smaller than for an insulated 
surface at the same Mach number; and by equation (71), 
the effect of withdrawing heat from the fluid 
through the solid surface is to increase Rs and to stab- 


ilize the flow, as compared with the flow over an insulated 
surface at the same Mach number (fig. 6). When the veloc- 
ity of the free stream at the “edge” of the boundary layer 
is supersonic, the laminar boundary layer is . completely 
stabilized if the rate at which heat is withdrawn through 
the solid surface reaches or exceeds a critical value that 
depends only on the Mach number, the Reynolds number, 
and the properties of the gas. The critical rate of heat 

transfer is that for which the quantity J- (^°f) * s - su (S - 

ciently large negatively near the surface (see equation (74)) 


so that (1— 2X)r(c) = 0.580 when c=c Q =l — ^(sections 4 

ui 0 

and 5). Although detailed stability calculations for super- 
sonic flow over a noninsulated surface have not been carried 
out, the function (1 — 2X)r(c) has been computed for non- 
insulated surfaces at fl/ 0 = 1.30, 1.50, 2.00, 3.00, and 5.00 by 
a rapid approximate method (appendix C) . The correspond- 
ing estimated values of Rf ■were calculated from equa- 

min 

tion (71), and in figure 7 these values are plotted against T u 
the ratio of surface temperature (deg abs.) to free-stream tem- 
perature (deg abs.). At any given (Mach number greater 

than unity, the value of Re increases rapidly as c 0 — U — iji 

07 miit aVaq 


when c 0 differs only slightlv from 1 — -tt-.* the stability of 

the laminar boundary is extremely sensitive to thermal 
conditions at the solid surface. At each value of fl/ 0 (> 1 , 
there is a critical value of the temperature ratio T x for which 
Rs r — > ®> . If T x ii T x . the laminar boundarv laver is 

c ‘ min. cr 

stable at all Reynolds numbers. The difference between the 
stagnation-temperature ratio and the critical-surface- 
temperature ratio, which is related to the heat-transfer 
coefficient, is plotted against Mach number in figure 8. 
Under free-flight conditions, for Mach numbers greater 
than some critical Mach number that depends largely on 
the altitude, the value of T s — T x ^ is within the order of 
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magnitude of the difference between stagnation temperature 
and surface temperature that actually exists because of heat 
radiation from the surface (references 14 and 15). In other 
words, the critical rate of heat withdrawal from the fluid 
for laminar stability is within the order of magnitude of the 
calculated rate of heat conduction through the solid surface 
which balances the heat radiated from the surface under 
equilibrium conditions. The calculations in appendix E 
show that this critical Mach number is approximately 3 at 
50,000 feet altitude and approximately 2 at 100,000 feet 
altitude. Thus, for Mo >3 (approx.) at 50,000 feet altitude 
and M 0 >2 (approx.) at 100,000 feet altitude, the laminar 
boundary-layer flow for thermal equilibrium is completely 
stable in the absence of an adverse pressure gradient in the 
free stream. 

When there is actually no heat conduction, through 
the solid surface, the limit of stability of the laminar 
boundary layer .depends only on the free-stream Mach 
number, that is, on the extent of the "aerodynamic- heating’' 


of the order of 


Mi 


» ('bu* 

\dy*; 


near the solid surface. A good 


indication of the influence of the free-stream Mach number 
on the distribution of p ^ across the boundary layer for an 
insulated surface is obtained from a rough estimate of the 
location of the point at which ^ (p reaches a positive 

maximum (or J-jp (p vanishes^. Differentiating the 

dynamic equations of mean motion twice and making use of 
the energy and continuity equations yields the following- 
result for an insulated surface: 


[d 3 fw'Y] _ 6 2 «) 2 
Jf\TJ J, 2 Ti m+2 


(76) 


where b = 5 


lie’ 


From equations (75) and (76) the value 


\v^x* 

of c at which (^t) van i s ^ es ) or ^ ^\T) reac ^ es 
mum, is given roughly for air by 


a maxi- 


, r & 

(w'X 

\J¥ 

W_l } 

~d 3 i 

'w'X 

df v 

,2yJ, 


7) 2 


(1 + 0 . 202 §Mq) 1 A ‘ 


(77) 


in which w/ » ■ ■ ) (appendix B). In other words, the 

point in the fluid at which (^t) stains a maximum moves 

farther out from the surface as. the Mach number is in- 
creased — at least in the range 0^M 0 ^4.5 (approx.); there- 
fore the value of c for which ^ vanishes and the value 
of c— c 0 for which (1— 2\)v(c) reaches the value 0.580 both 


increase with the Mach number (fig. 3 (a)). By equation 
(71), the value of R»a m(n for the laminar boundary-layer flow 

over an insulated surface decreases as the Mach number 
increases and the flow is destabilized, as compared with the 
Blasius flow (fig. 5). 

C. RESULTS OF DETAILED STABILITY CALCULATIONS FOR INSULATED AND 
NONINSULATED SURFACES 

From .the results of the detailed stability calculations for 
several representative cases (figs. 4 to 6), a quantitative esti- 
mate can be made of the effect of free-stream Mach number 
and thermal conditions at the solid surface on the stability of 
laminar boundary-layer flow. For the insulated surface, the 
value of Re is 92 when M 0 = 1.30 as compared with a value 

cr min 

of 150 for the Blasius flow. For the noninsulated surface at 
Mo=0.70, the. value of Re is 63 when 7) = 1.25 (heat 

added to fluid), Re. r =126 when 7\=1.10 (insulated sur- 

cr min 

face), and Re =5150 when 7\=0.70 (heat withdrawn 

cr min 

from fluid). Since R x m « 2.257?s 2 (the value of 6 \ f i=r~jwhich 

V v a *x* 

is proportional to the skin-friction coefficient, differs only 
slightly from the Blasius value of 0.6667), the effect of the 
thermal conditions, at the solid surface on R / is even more 
pronounced. The value of R/ is 60 X 10 s when 7\=0.70 and 
Mo=0.70, as compared with a value of 51 X 10 s for the 
Blasius flow .(T^l and M 0 = 0). For the insulated surface 
the value of R x * declines from the Blasius value for 

"min 

M 0 =0 to a value of 19X10 3 at M 0 =1.30. The extreme 
sensitivity of the limit of stability of the laminar boundary 
layer to thermal conditions at the solid surface when Tf-O 
is accounted for by. the fact that c 0 is small when TXl and 
M 0 <1 (or Mo is not much greater than unity) and 

Re (equation (71)). Small changes in c 0 , therefore, 

"min c 0 

produce, large changes, in 7?« . In addition, when T’i’O, 

CT min 

small changes in the thermal conditions at the solid surface 
produce appreciable changes in (yf) (equation (74)) and, 
therefore, in the value of c 0 . 

Not only is the value of I?« affected by the thermal 

/ cr min 

conditions at the solid surface and by the free-stream Mach 
number but the entire curve of a$ against 7? s for the neutral 
disturbances is also affected. (See figs 4 (k) and 4 (1).) 
When the surface is insulated (and M 0 XO), or heat is added 
to the fluid (7 1 i =1.25), as 7? 9 — > <» along the upper 

branch of. the curve of neutral stability. In other words, 
there is a finite range of unstable wave lengths even in the 
limiting case of an infinite Reynolds number (inviscid fluid). 
However, a^O as Ru— »«> for the Blasius flow, or when heat 
is withdrawn from the fluid. This behavior is in complete 
agreement with the results obtained in section 4 and. in 
reference. 8. 
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A comparison between the curves of at against Re for 
7,= 1.25 and 71 = 0.70 at fl/ 0 =0.70 shows that withdrawing 
heat from the fluid not only stabilizes the flow by increasing 
Re but also greatly reduces the range of unstable wave 

€T min ^ 

numbers (at). On the other hand, the addition of heat to 
the fluid through the solid surface greatly increases the range 
of unstable wave numbers. 

It should also be noted that for given values of at, c, 
and Re the time frequencies of the boundary-layer disturb- 
ances in the high-speed flow of a gas are considerably greater 
than the frequencies of the familiar Tollmien waves observed 
in low-speed flow. The actual time frequency «* expressed 
nondimensionally is as follows: 

ll*vn* _ Cat 
(w 0 *) 2 2xi?s 

For given values of c, a s , and Rt the frequency increases as 
the square of the free-stream velocity. 

D. INSTABILITY OF LAMINAR BOUNDARY LAYER AND TRANSITION TO 
TURBULENT FLOW 

The value of Re obtained from the stability analysis 

cr min 

for a given laminar boundary-layer flow is the value of the 
Reynolds number at which self-excited disturbances first 
appear in the boundary layer. As Prandtl (reference 12) 
carefully pointed out, these initial disturbances are not 
turbulence, in any sense, but slowly growing oscillations. 
The value of the Reynolds number at which boundary- 
layer disturbances propagated along the surface will be 
amplified to a sufficient extent to cause turbulence must be 
larger than Re in any case; for the insulated flat-plate 

flow at low speeds and with no pressure gradient, the tran- 
sition Reynolds number Rt tr is found to be three to seven 
times as large as the value of i?» (references 6 and 7). 
The value of Re depends not only on Re but also on 

the initial magnitude of the disturbances with the most 
‘‘dangerous” frequencies (those with greatest amplification), 
on the rate of amplification of these disturbances, aDd on the 
physical process (as yet unknown) by which the quasi- 
stationary laminar flow is finally destroyed by the amplified 
oscillations. (See, for example, references 16 and 17.) The 
results of the stability analysis nevertheless permit certain 
general statements to be made concerning the effect of 
free-stream Mach number and thermal conditions at the 
solid surface on transition. The basis for these statements is 
summarized as follows: 

(1) In many problems of technical interest in aeronautics 
the level of free-stream turbulence (magnitude of initial 
disturbances) is sufficiently low so that the origin of transition 
is always to be found in the instability of the laminar bound- 
ary layer. In other words, the value of Re crmtn I s absolute 
lower limit for transition. 

(2) The effect of the free-stream Mach number and the 
thermal conditions at the solid surface on the stability limit 


(Rtcr^J) is overwhelming. For example, for M a =0.70, the 
value of Re„ when 71 = 0.70 (heat withdrawn from fluid) 

cr min. 

is more than 80 times as great as the value of Re„ mill when 


7\=1.25 (heat added to fluid). 

(3) The maximum rate of amplification of the self-excited 
boundary-layer distur bances propagated along the surface 
varies roughly as 1 (This approximation agrees 

closely with the numerical results obtained by Pretsch 
(reference 18) for the case of an incompressible fluid.) The 
effect of withdrawing heat from the fluid, for example, is 
not only to increase Re,. . and stabilize the flow in that 

*- m ut 

manner but also to decrease the initial rate of amplification 
of the unstable disturbances. In other words, for a given 
level of free-stream turbulence, the interval between the 
first appearance of self-excited disturbances and the onset of 
transition is expected to be much longer for a relatively 
stable flow, for which is large, than for a relatively 


unstable flow, for which Re C t nbl is small and the initial rate 
of amplification is large. 

On the basis of these observations, transition is delayed 
(Re lr increased) by withdrawing heat from the fluid through 
the solid surface and is advanced by adding heat to the fluid 
through the solid surface, as compared with the insulated 
surface at the same Mach' number. For the insulated surface, 
transition occurs earlier as the Mach number is increased, as 
compared with the flat-plate flow at very low Mach numbers. 
\,hen the free-stream velocity at the edge of the boundary 
layer is supersonic, transition never occurs if the rate of heat 
withdrawal from the fluid through the solid surface reaches^ 
or exceeds a critical value that depends only on the Mach 
number (section 6B and figs. 7 and S). 

A comparison between the results of the present analysis 
and measurements of transition is possible only when the 
free-stream pressure gradient is zero or is held fixed while the 
f re e-stream Maeh number or the thermal conditions at the 
solid surface are varied. Liepmann and Fila (reference 19) 
have measured the movement of the transition point on a 
flat plate at a very low free-stream velocity when heat is. 
applied to the surface. They found by means of the hot-wire 
anemometer that declined from 5X10* for the insulated 
surface to a value of approximately 2X10 S for 71=1.36 when 


the level of free-stream turbulence \l 


/( w * 7 ) 2 




was 0.17 percent, 


or to a value of 3X10 S when 


V ( 


, ===0.05 percent and 

Uo*j 2 


71=1.40. The value of R Slr declines from 470 (approx.) to 
300 (approx.) in the first case and to 365 in the second. 

Frick and McCullough (reference 20} observed the varia- 
tion in the transition Reynolds number when heat is applied 
to the.upper surface of an NACA 65, 2—016 airfoil at the nose 
section alone, at the section just ahead of the minimum pres- 
sure station, and for the entire laminar run. When heat is 
applied only to the nose section, the transition Reynolds 
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number (determined by total-pressure-tube measurements) 
was practically unchanged. Near the nose, :Rj e 
and the strong favorable pressure gradient in the region of 
the stagnation point stabilizes the laminar boundary layer 
to such an extent that the addition of heat to the fluid has 
only a negligible effect. When heat is applied, however, to 
the section just ahead of the minimum pressure point, where 
the pressure gradients are moderate, the transition Reynolds 
number Re tr declined to a value of 1190 for T x ^ 1.14, com- 
pared with a value of 1600 for the insulated surface. When 
heat is applied to the entire laminar run, R Blr declined to a 
value of 1070 for 7)*= 1.14. 

It would be interesting to investigate experimentally the 
stabilizing effect of a withdrawal of heat from the fluid at 
supersonic velocities. At . any rate, on the basis of the 
results obtained in the experimental investigations . of - the 
effect of heating on transition at low speeds, the results of 
the stability analysis give the proper direction of this effect. 


7. STABILITY OF LAMINAR BOUNDARY-LAYER FLOW OF 
A GAS WITH A PRESSURE GRADIENT IN THE DIRECTION 
OF THE FREE STREAM 


For the case of an incompressible fluid, Pretsch (reference 9) 
has shown that even with a pressure gradient in the 
direction of the free stream, the local mean-velocity distri- 
bution across the boundary laj'er completely determines the 
stability characteristics of the local laminar boundary-layer 
flow at large Reynolds numbers. From physical considera- 
tions this statement should apply also to the compressible 
fluid, provided only the stability of the flow in the boundary 
layer is considered and not the possible interaction of the 
boundary layer and the main “external” flow. Further 
study is required, to settle this question. 

If only the local mean velocity-temperature distribution 
across the boundary layer is found to be significant for 
laminar stability in a compressible fluid, the criterions ob- 
tained hi the present paper and in reference 8 are then im- 
mediately applicable to laminar boundary-layer gas flows in 
which there is a free-stream pressure gradient. The quanti- 
tative effect of a pressure gradient on laminar stability could 
be readily determined by means of the approximate estimate 
of Rscr mi „ (equation (70)) in terms of the distribution of 


the quantity p ^ across the boundary layer. 


Such calcula- 


tions (unpublished) have already been carried out by. 
Dr. C. C. Lin of Brown University for the incompressible fluid 
by means of the approximate estimate of given in 

reference 5, part III. 

In any event, the qualitative effect of a free-stream pres- 
sure gradient on the local distribution of p across the 

boundary layer is evidently the same in a compressible fluid 
as in an incompressible fluid. If the effect of the local 
pressure gradient alone is considered, the velocity distribu- 
tion across the boundary layer is “fuller” or more convex for 
accelerated than for uniform flow and, conversely, less 


convex for decelerated flow. Thus, from the results of the 
present paper the effect of a negative pressure gradient on 
the laminar boundary-layer flow of gas is stabilizing, so far 
as the local mean velocity-temperature distribution is con- 
cerned, while a positive pressure gradient is destabilizing. 
For the incompressible fluid, this fact is well-established by 
the Rayleigh-Tollmien criterion (reference 3), the work of 
Heisenberg (reference 1) and Lin (reference 5), and a mass of 
detailed calculations of stability limits from the curves of a 
against R for the neutral disturbances. These calculations 
were recently carried out by several German investigators 
for a comprehensive series of pressure gradient profiles. 
(See, for example, references 9 and 21.) 

Some'idea of the relative influence on laminar stability of 
the thermal conditions at the solid surface and the free-stream 
pressure gradient is obtained from the equations of mean 
motion. At the surface, 



— Po 


■ — * du o* 

IP 


(78) 


or 


KC 


dw\~ 
p dy)_ 


m4- 1 


2V 


7Vu>.'- 


1 5 2 du a * 

dx* 


(79) 


J-AJ. >_> i 

l\ S 2 du 0 * \ 

lt*W s2 - 8a3 7 


the distribution of p j- is sensitive to the 


In a region of small or moderate pressure gradients 

dw 
dy 

thermal conditions at the solid surface. For example, the 
chordwise position of the point of instability of the laminar 
boundary layer on an airfoil with a flat pressure distribution 
is expected to be strongly influenced by heat conduction 
through the surface. (See reference. 20.) For the insulated 
surface, the equations of mean motion yield the following 
relation (appendix D), which does not involve the pressure 
gradient explicitly: 



= v(m-j- 1) (y — 1) il/ 0 2 


WY 
l V ' 


>o 


(80) 


The effect of “aerodynamic heating” at the surface opposes 
the effect of a favorable pressure gradient so far as the dis- 
tribution of p across the boundary layer is concerned 

(equations (79) and (80)). The relative quantitative influ- 
ence of these, two effects on laminar stability can only be 
settled by actual calculations of the laminar boundary-layer 
flow in a compressible fluid with a free-stream pressure 
gradient. A method for the calculation of such flows over 
an insulated surface is given in reference 22. 

When the local free-stream velocity at the edge of the 
boundary layer is supersonic, a negative pressure gradient 
can have a decisive effect on laminar stability. The local 
laminar boundary-layer flow over an insulated surface, for 
example, is expected to be completely stable when the mag- 
nitude of the local negative pressure gradient reaches or 
exceeds a critical value that depends only on the local Mach 
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number and the properties of the gas. The critical magni- 
tude of the pressure gradient is that which makes the quantity 


A 

dy 


{ dw 
\ dy, 


sufficiently large negatively near the surface so that 


-[ 1 - 


•2 


T 2 

(v/f 


d fyf 
dy\T 


=0.580 


when c=l — -yy* 

It has already been shown in the present paper than when 
il/o)>3 (approx.) the laminar boundary-layer flow with a 
uniform free-stream velocity is completely stable under free- 
flight conditions when the solid surface is in thermal equi- 
librium, that is, when the heat conducted from the fluid to 
the surface balances the heat radiated from the surface (sec- 
tion 6B). The laminar boundary-layer flow for thermal 
equilibrium should be completely stable for say, 

where M s <d . 3 if there is a negative pressure gradient in the 
direction of the free stream. Favorable pressure gradients 
exist over the forward part of sharp-nosed airfoils and bodies 
of revolution moving at supersonic velocities, and the limits 
of stability (Re cr ) of the laminar boundary layer should 
be calculated in such cases. 


CONCLUSIONS 


From a study of the stability of the laminar boundary 
layer in a compressible fluid, the following conclusions were 
reached : 

1. In the compressible fluid as in the incompressible fluid, 
the influence of viscosity on the laminar boundary-layer flow 
of a gas is destabilizing at very large Reynolds numbers. 
If the free-stream velocity is subsonic, any laminar boundary- 
layer flow of gas is unstable at sufficiently high Reynolds 
numbers. 

2. Regardless of the free-stream Mach number, if the 
product of the mean density and the mean vorticity has an 


extremum 


d ( dw\ 
dy v dy) 


vanishes^ for some value of w>l — 

(where w is the ratio of mean velocity component parallel 
to the surface to the free-stream velocity, and where ilia is 
the free-stream Mach number) the flow is unstable at suffi- 
ciently high Reynolds numbers. 

3. The actual limit of stability of laminar boundary-layer 
flow, or the minimum critical Reynolds number Rt„ miri , is 
determined largely by the distribution of the product of the 
mean density and the mean vorticity across the boundary 
layer. An approximate estimate of Rt CTmin is obtained that 
serves as a criterion for the influence of free-stream Mach 
number and thermal conditions at the solid surface on 
laminar stability. For zero pressure gradient, this estimate 
reads as follows: 


R> 


[W 6 


Ti 


Cq 


vR 


c 0 ) 2 


where T is the ratio of temperature at a point within the 
boundary layer to free-stream temperature, T 1 is the ratio 
of temperature at the solid surface to the free-stream tem- 


perature, and c 0 is the value of c (the ratio of phase ve- 
locity of disturbance to the free-stream velocity) for which 
(1— 2X)r=0.580. The functions c(c) and X(c) are defined as 

follows : 



T 2 b 
/ dw \ 3 Si? 

Ad n ) 


1 _ dw\ 
Tdy) 



1 


where 


V 


nondimensional distance from surface 



4. On the basis of the stability criterion in conclusion 3 
and a study of the equations of mean motion, the effect of 
adding heat to the fluid through the solid surface is to 
reduce R$a- mill and 1° destabilize the flow, as compared with 
the flow over an insulated surface at the same XIach number. 
Withdrawing heat through the solid surface has exactly the 
opposite effect. The value of i? fcrmfa for the laminar boundary- 
layer flow over . an insulated surface decreases as the 
Mach number increases, and the flow is destabilized, as 
compared with the Blasius flow at low speeds. 

5. When the free-stream velocity is supersonic, the lami- 
nar boundary layer is completely stabilized if the rate at 
which heat is withdrawn from the fluid through the solid 
surface reaches or exceeds a certain critical value. The 
critical rate of heat transfer, for which Rhr mta -^ <n , is 

d, f dw\ 

that which makes the quantity sufficiently large 

negatively near the surface so that [1— 2X(c)J v(c) =0.580 

when c=c Q =l — p- Calculations for several supersonic 
iV-lQ 

XIach numbers between 1.30 and 5.00 show that for AJ' 0 >3 
(approx.) the critical rate of heat withdrawal for laminar 
stability is wi thin the order of magnitude of the calculated 
rate of heat conduction through the solid surface that bal- 
ances the heat radiated from the surface under free-flight 
conditions. Thus, for X/ 0 > 3 (approx.) the laminar boundary- 
layer flow for thermal equilibrium is completely stable 
at all Reynolds numbers in the absence of a positive (adverse) 
pressure gradient in the direction of the free stream. 

6. Detailed calculations of the curves of wave number 
(inverse wave length) against Reynolds number for the 
neutral boundary-layer disturbances for 10 representative 
cases of ins ulated and noninsulated surfaces show that also 
at subsonic speeds the quantitative effect on stability of 
the thermal conditions at the solid surface is very large. 
For example, at a XIach number of 0.70, the value of Rt C r m!x 
is 63 when 27 = 1.25 (heat added to fluid), Recr mix = 126 
when 27 = 1.10 (insulated surface), and Rg„ mu = 5150 when 
37= 0. 70 (heat withdrawn from fluid). Since R X .^2.25R } 2 , 
the effect on R x -cr miyi is even greater. 
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7. The results of the analysis of the stability of laminar 
boundary-layer flow by the linearized method of small per- 
turbations must be applied with care to predictions of 
transition, which is a nonlinear phenomenon of a different 
order. Withdrawing heat from the fluid through the solid 
surface, however, not only increases but decreases 

the initial rate of amplification of the self-excited disturb- 
ances, which is roughly proportional to 1 l-\Rhr mi J addition 
of heat to the fluid through the solid surface has the opposite 
effect. Thus, it can be concluded that (a) transition is 
delayed (Re lr increased) by withdrawing heat from the fluid 
and advanced by adding heat to the fluid through the solid 
surface, as compared with the insulated surface at the same 
Mach number, (b) for the insulated surface, transition 
occurs earlier as the Mach number is increased, and (c) when 
the free-stream velocity is supersonic, transition never 
occurs if the rate of heat withdrawal from the fluid through 
the solid surface reaches or exceeds the critical value for 
which Rhr min ~* “ ■ (See conclusion 5.) 

Unlike laminar instability, transition to turbulent flow in 
the boundary layer is not a purely local phenomenon but 
depends on the previous history of the flow. The quantita- 
tive effect- of thermal conditions at the solid surface on 
transition depends on the existing pressure gradient in the 
direction of the. free stream, on the part of the solid surface 
to which heat is applied, and so forth, as well as on the 
initial magnitude of the disturbances (level of free-stream 
turbulence) . 

A comparison between conclusion 7 (a), based on the 
results of the stability analysis, and experimental investi- 


gations of the effect of surface heating on transition at low 
speeds shows that the results of the present paper give the 
proper direction of this effect. 

8. The results of the present study of laminar stability 
can be extended to include laminar boundary-layer flows 
of a gas in which there is a pressure gradient in the direction 
of the free stream. Although further study is required, it 
is presumed that only the local mean velocity-tempera ( ure 
distribution determines tl^e stability of the local boundary- 
layer flow._ If that should be the case, the effect of a i>ressure 
gradient on laminar stability 1- could be easily calculated 
through its effect on the local distribution of the product 
of mean density and mean vorticity across the boundary 
la} r er. 

When the free-stream velocity at the “edge” of the bound- 
ary layer - is supersonic, by analogy with the stabilizing 
effect of a withdrawal of heat from the fluid, it is expected 
that the laminar boundary-layer flow is completely stable 
at all Reynolds numbers when the negative (favorable) 
pressure gradient reaches or exceeds a certain critical value 
that depends only on -the Mach number and the properties 
of the gas. The laminar boundary-layer flow over a surface 
in thermal equilibrium should be completely stable for 
M 0 >A/„ say, where if there is a negative pressure 

gradient in the direction of the free stream. 


Langley Memorial Aeronautical Laboratory, 
National Advisory Committee for Aeronautics, 
Langley Field, Va., September 5, 1940. 



APPENDIX A 


CALCULATION OF INTEGRALS APPEARING IN THE INVISCID SOLUTIONS 


In order to calculate the limits of stability of the laminar boundary layer from relations ( 21 ) to (29) between the values 
of phase velocity, wave number, and Reynolds number, it is first necessary to calculate the values of the integrals K u 
Hi, Hi, Ah, il/ 3 , A r 3 , and so forth, which appear in the expressions for the inviscid solutions <?i(y) and y> 2 (y) and their 
derivatives at the edge of the boundary layer. These integrals are as follows (equations (13), (9), and (10)): 


Hi (c) = 
Kdc) = 


'v'-iii'-cy 


dy 


ei T—M 0 2 (w—c) 2 ] 
i — -dy 


(w—c) 2 


N 2 (c) — KiHi— K,= 


,fl T—M 0 2 (w—cY 


, C l iU'—C ) 2 , T! , 

2 dy — 77 . dy = H 2 (c) 

Jft {W—C) Jifi J- 


M,(c)=H 2 H l -H 3 = ( 

tttt tz f" T—MYiw—cY 


S! (w— c ) 2 , C^-T— MY(v}— c'Y , r* (w—c) 2 


{ w—c Y 


l v 

dy 

J S'. 


-dy 




and so forth. 

Terms of higher order than a 3 in the series expressions for er and are neglected. When a<(l, the error involved is 

. ft 

small because the terms in the series decline like Even for «)> 1 , however, this approximation is justified, at least 

for the values of c that appear in the stability calculations for the 10 representative cases selected in the present paper. 

I r c 3 ~\ k ~ 1 

For example, the leading term in R. P. AVhM. where 1—2,3, . . ., is approximately ^ multiplied by 


the leading term in R. P. N 3 (c). The quantity in the brackets is at most 0.12 in the present calculations; for example, 
R. P. N s (c) ~ 0.06 R. P. A r 3 (e). Moreover, R. P. A r 2t (c)«(l — c ) R. P. A T 2 i + 1 '(c). Similar approximate relations exist 

between R. P. il I 2k (c) and R. P. M 3 (c); and, in addition, R. P. M 3 (c)~( 1 — c) ^ R. P. A r 3 (c) —0.015 R. P. N 3 (c), at most. 


The only integral for which the imaginary part is calculated is K\(c). At the end of this appendix, it is shown that 
the contributions of the imaginary parts of f? 2 , M it and A T 3 are negligible in comparison with the contribution of I. P. Ki(c). 


GENERAL FLAN OF CALCULATION 

The method of calculation adopted must take into account 

the fact that the value of ^(p Hy) at *-be P 0 ^ y=Vc 

where w=c, strongly influences the stability of the laminar 
boundary layer. Accordingly, the integrals are broken into 
two parts; for example, 


Ki 


(■Vi T fj; T 

i(c)= J, ■&=& dy+ L 

=A II (c)+R 12 (c)-A/ 0 2 

where yH>Vc- The integral Kn(c), which involves 
\jHj C P ~3dy) ’ * s ca ^ cu ^ a t e< i ver T accurately, whereas Kn(c) 

is calculated by a more approximate method as follows: 


r y - t 

° ! ' )= J„ <a=& iv 


(Al) 


This integral is evaluated as a power series in c. The velocity 
profile w(y) is approximated by a parabolic arc plus a 
straight-line segment for purposes of integration. In the more 
complex integrals H 2 , and A r 3 , the indefinite integrals 
r? x Ei. T 

i-s dy and 7 - 2 dy are evaluated by 21 or 41 

J yi (w-cf J ( w—c ) 2 J 

point numerical integration by means of Simpson’s rule. 
The values of w(y ) are read from the velocity profiles of 
figures 1 and 2 . The value of y—y^—a is 0.40 in the present 
series of calculations; this value is chosen so that the point 
y=Vi is never too close to the singularity at y=yc Take 



T J 
(w—c) 2 y 


(A 2) 


T 


dy 


The integral A* u (c), or the indefinite integral j* ^ w _ c ^ 

that appears in H 2 , il/ 3 , and N s , is evaluated by expanding 
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the integrand in a Taylor’s series in y—y c and then integrat- 
ing the series term by term. The path of integration must 
be taken below the point y=y c in the complex y-plane. 

Instead of calculating the values of the velocity and tem- 
perature derivatives w c in) and T c lK) directly, it is simpler 
to relate these derivatives to their values at the surface by 
Taylor’s series of the form 


w c in) = Wi w + Wi <K+I: 1 (y e — Vi) -f 


Wl 


( n+2 ) 


21 


i.Vc-yi) 2J r 


The derivatives., at the surface w>i <B> and Ti U) are calculated 
from the equations of mean motion (appendix B) . 

The integral K u (c), for example, is finally obtained as a 
power series in y c —y \=<r and in y } —y c =a — a, plus terms 
involving log <r. The phase velocity c is related to <r by 




A s 


3 = Wj' ^<T + ™ rS " 


where 


Ai 


Wt 


00 




Terms up to the order of a 1 are retained in order to include 
all terms involving w 1 v11 . 


DETAILED CALCULATIONS 


In order to illustrate the method , the evaluation of K x ( c ) 
is given in some detail, as follows: 

(1) Evaluation of K\{c): 


K\ 


(a) Define 
Now 


fvt J’ 
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T T 

(w— c) 2 (w/) 2 (y-y c )V 


where 


'P(y)=i+^P (y-yc)+j(^p-(y-yc) 2 + 


T . 

The function is developed in a Taylor’s series around the 
point w=c as follows: 

.+(£), {v ~ v ‘ )+ 

where 
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where 

Vi— yc=\yi—yde~ iT 

Vi—Vc= (: Vf—yi ) — (Vc—Vi) =a-cr 

°=y—yi 

The coefficients are expressed in terms of derivatives 

of T and w at y=y x as follows: 

Define 


Then 
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=Jk{yi)+f/{yi)(y c -yi)+— 2 Y^ {y-yiYA. ■ ■ 

(The method adopted for the calculation of /* n (y l ) from the 
velocity and temperature derivatives w x u) and T X U) is given 
at the end of this appendix.) 

From the expression for if n (c), 

I. P. K n (c)=I. P. KM 

=vrfi (y c ) 


/i(2/i) + tr /i / (2/i)+ • ■ • /i v (i/i)J 


and 


T x 
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where 

<r=y—yi 


ft— ~dk+ifa(yi) O^Jc 


c ft =s*+-wrwTvf-^*+Do(yi) 0^/r<5 (s 5 =0) 

Sz=aj 2 {y l )+a%(y l )+a%{y 1 )+a%{y l )+a%{y l ) + • . . 
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Si = aji (yi) J ra% f (y 1 )+ a 3 j '/ (y : ) +a%'(y 1 )+ . 
[2afs(.Vi) + 3a% (y t ) + 4g 3 / 5 (i/ ! ) + 5a 4 /a(l/i) + 


■ 1 


laft'iyi) +a*f,"(y l )+a?j i "(y l ) + . . . ]- 

[2af»{yi) + 3a%'(yi)+4a?f s , (y 1 ) + . . . ] + 
[2>af i (y 1 ) + Qa%(y l ) + 10a%{y 1 )+ . . . ] 

s i =\[aj, n, {yi) + a%"’{y 1 )+ ...]- 

\[2aj/'{y l )+3a%"(y 1 )+ ...]+ 

[3af/(y l )-{-e,a%'(y 1 )+ ...J— 

[4q/ 5 (yi) + 10a 2 / 6 (i/i)-i- . . . ] 

[«/**% i)+ • - - }-ll2af 3 '"(yi)+ . . . ] + 
\[3aji"{yi)+ . . . ]-[4a/ 5 '(y 1 )+ ... 1+ 

[5a/ s (t/i) + . . . ] 

<£s = — X) (r+1)! ^ l - T (^o= 1-0) 

A t =^r (o=0.40) 

(b) Define 

pi.o j’ 

J 0.40 (W — c) 2 d(.y — Vl) 

=X a k (.k+ l)c l 

k = 0 

where 

m.o ]y 

a ‘~ 


The velocity profile w(y) is approximated by a parabolic 
arc in the interval 0.40 ^y— yi^y 3 — y x and by a straight 
line (tt>=Constanfc=i 0 (jf*))in the interval y 3 — y t 25y— y^ 1.0. 
The value of y 3 is determined by imposing the condition 
that the area under the parabolic-arc straight-line 
segment equals the area under the actual velocity profile 
w(y) in the interval 0.40^y— zqSs 1.0. The parabolic 
arc w=l+m(y—y 1 )+n(y—y i y is determined by the 
following conditions: when y=y 4 <0, 

w=l 


w '=0 


Therefore 

&k = T\{Ik+'i J rJk+‘l) ^(2"l — 1) “ ~2 + 

M 0 2 (I t +J t ) 

where 

r = f*-" d(y—yi) 

4 Jo. 40 V? 

and 

T _ f L0 d(y— yO_l — fa— yi) 

[W(jft)l* My*)]* 

I t is evaluated by approximating w(y) by a parabolic arc 
as follows: 




7 X - J- in f VZ— m— 2n(y— yQ l* 

-\ A ]_-\ r A+m+2n(,y— yOJo.M 

1 r ra-\-2n{y—y-i) "j 1,1-1 

(k— i)A \_[l+ m (y— yi).4-n(y— yiYY' 1 !^ 

2k— Z 4(— n) r 

2k— 2 A lb ~ l - - 


where A—m-—41n. 

As a control in the calculation of the series expression 

co 

yi a k (k+l)c k for Ku(c), use is made of the fact that, from 

fc =0 

the definition of I t and J t) 


lim (J*+ J*) = — 
* — ‘ ■“ k 


and therefore 


!5|r] 


lim 


.Cw> 


w (Vj) k-\-l 


The remainder after N terms in the series for K u {c) is given 
approximately by 

f(A r -[-l) term] 



The real part of Ki(c) is obtained by combining the results 
of (a) and (b); that is, 

R. P. K, (c)=R. P. [z„ (c) +^;-]+Z I2 (c) -ff/ 0 2 
(2) Evaluation of Hi(c): 


II 


(c)= P 
J Hi 


{w—c)' 


dy 


when y=yi and y } — yi=0.40, 

w=w{yj) 1 

where w(y } ) is read off the velocity profile of figures 1 and 2. 
The value of y 4 is chosen so that the parabolic arc fits the 
velocity curve w{y) closely over the widest possible range. 
For cr= 1 , 

T= 2\— [(Ti— 1) av] M/w 2 


The integrand of this integral is free of singularities in the 
region of the complex y-plane bounded by y=y\ and y=y 2 ; 
therefore, Hi{c) is evaluated by purely numerical integration. 
The actual procedure employed for the calculation of in- 
tegrals of this type is as follows: (The integral Hi(e) serves 
as an illustration.) 

(a) Define 

Si{c)=^^ pW 2 dy— 2c J' plffdrj -j- C 2 J ' pd’qj 
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where 


and 


■n =y* 


V* 

=iv'J5- 

V V 0 *x* 


(b) With the approximation that the viscosity varies linearly 
with the absolute temperature, the velocity to is the same 
function of the nondimensional .stream function f as in the 
Blasius flow; that is, 

where { is defined by the relation dt—pwdri (appendix B). 
From these relations 

P w n dr\=[w($)) n - 1 d$=[w B (r lB )] n d?i B 
since d{=w B dv B - Moreover, 


d r 


where 


dr, =MT) =TiWB)d,lB 

' ] w 


T{W B ) = T, - [ ( T x - 1 ) - : ^ Mr 2 


7—1 

cy ArlQ UJb 2 

for <r= 1 . 

(c) Finally, from the relations given in (b), 


AIq 2 Wb 2 


w B 2 dr] B — 2c w B dijB+c 


u n 


where b Q is the value of 5 y for the Blasius flow. For 


I v a *x* 

the insulated surfaces, b 0 , which is somewhat arbitrary, was 
chosen as 5.60; whereas for the noninsulated surfaces, 
6 0 =6.00. (The value of w B at ij s =5.60 is 0.9950; when 
i7 S =6.00, " 103 = 0.9975. The value of b for the insulated 
surfaces is the value of ij at which w=0.9950 ; whereas b for 
the noninsulated surfaces is the value of _ij for which 
w=0.9975.) The advantage of this procedure is that the 

w B n dr} B are calculated once and for all and the 
o 

value of Hi(c) depends only upon the values of b and c. 
In fact, 


since 


and 


Also, 


HM 


730 


r 

Jo 


} =0.6667 
v 0 *x*J B 

w B dy B =b a — 1.730 


and 


C b r b ‘ 
b= di ;= T drts 

Jo Jo 


= 6 0 +1.73(T 1 -l) + 0.6667 ^ M 0 l 
=6 0 + 1.73 j^(Ti— 1) ~"^2 — fl^o 2 J+2. 


3967 -'- 2 1 AfJ 


See appendix B. (Incidentally, the last relation shows the 
effect of free-stream Mach number and thermal conditions 
at the solid surface on the “thickness” of the boundary 
layer.) 

(3) Evaluation of H 2 (c): 


H, 


«-x 


T—M 0 2 (w—c ) 2 


m (w—cY 

Vi J 1 


, C y (w—cY 2 J 

dy j —j’ dy 


tn 

C V1 t 7 c* {w-cY j 1r , enw-cY , , 

J, ~r~ dy dy 


Define 


vt 7 1 


, i‘ 1 ' (w-c)' 2 , 
,Aw--i? dy ] „ -I— * 


H h (c) 

. . . H 2i (c) = J] dy dy 

(a) The integral Hz 2 (c) is evaluated by methods similar to 
those already outlined for the evaluation of Hi (c ) . Thus 

H22{c)= £Sl^^- dydy 

J 'vt ru rut rv rut rv 

dy pw 2 dy-2c\ dy pwdy+c 2 I dy p dy 
u i Jy i Jv i Jin Jvi Jin 


c h 

c 2 j o T ns dm 


W’s tfljs 


where 


T= 


= Tt -[( T x - 1) - Mo 2 ] w B -^l M 0 ‘w b ! 


The nine integrals in the expression for Ii 2l (c) are evaluated 
by numerical integration using Simpson’s rule. 

(b) Define. 

= £^ dv £,^ dy+ tiJtf iv L 


J I 9 (w—cY , 

„ {w-c) idy J„ T dy 


Define 






t Cv 

1 ^ Jin (W-C) 2 ^!,, 


Cm T Cv 

,M= , 2 dy dy 

12 jy,- (w-c) 2 T J 


(w—c)' 1 . 
j - dy 

(w— c) 2 
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The integral H 2 (c) is evaluated as follows: 


H, 


■ w = r isiWrv L ’ ii ' 

T , CHw—cYj 
T , f »*(«?— cl 2 . 


But 

and 
so that 


ryi T 

X, dy=K lt ( C ) 

-y~'dy=H<Cci 


r» ( W - C y 

J 


Define 


where 


and 


nu t rm ( w 

ry 2 r p 

plc)= L 0(y:c)iy 

i r 6 t 

= 5 * Jq .46 

f 6 v/, 0 f 6 w , L 2 f 6 dy 

~J n T dn 2c j, iT 77 1 C J, T 

r= 7 T 1 -[(r 1 -l)-^A/ 0 1 ] 


-cl 2 


T 


dy 


w-^MW 


The integral P(c) is evaluated by numerical integration using 
Simpson’s rule: the required values of w are read directly 
off the velocity profiles of figures 1 and 2 . Finally, 

mjc)=Ku(c) Hi(c)—P(c) 

The integral H 2n (c) is evaluated in exactly the same way as 
Ku(c) where 

I to— c) 2 


T 


and 


■-(w/Y{y-ycY (|T) 


^(t/)=i+^?(y~yc)+|^v(y-2/ c ) 2 + - - • 

R. P. [H 2 jc)] = (b (i a*-+c a az+d 0 a i +n 0 a s ) + 
tr( i 6 I u 2 -l-c I a 3 -i-d 1 a 4 —2a& 0 — 3a 2 Co—4a 3 do—5a 4n o)-r 
ff 2 (6 2 a 2 -|-C2a 3 — 2a5: — 3ct 2 Ci- — 4a 3 <fi — g-f- 3ctCo+ 6a 2 d 0 -r 10a s 7io) ~h 
c z (b ;i a 2 —2ab 2 —3a i C2-\-3aci-rQa t d 1 J rp 0 —4ad (I —10a s n 0 -\-a s )-r 
i7 i (—2ah z +3ac 2 J rp l — 4:ad 1 + g a J r5an 0 -\-a i ) J r 


w. 



i K 4 
a — cr (3 


l /o'Cyi) i i 

L 3 MvO 


A ' 4 \ 2 Tj\ 


+ 


UAM^MeIAM (a 

I 3 " IMVx)? ‘ ±MVi) \ 2 T x ) 3 Uy x ) \ ) 


where 
( it ’/) 2 


h 


6 


&i=y (mT 

6 l= (-ii±A) W) . 

‘•-(n 1 + ur) ta '>’-$£a 

„ My i) 

3 3/o(2/i) 

{ 1 Mh)M ( 2 /.) , 1 i 2 (.Vi) ( , 2 V \ 1 // (y, 

':)V 2 W S/ Q (y/ 


fli 


Cn = 


<V 


(3 L/a(y:)l 2 ‘ 4/0(2/:) 

1 /1 (2/1) 1 / < _Ti’\ 

‘9/0(2/!) 12 \ * rj 

(i //(yQ i fi(yi)fa'(yi ) , i 
(9/0(2 /:) 9 L/o(2/i)] 2 ‘12 


O'-f+OT 


, f 1 /1 (yi) r/o / (y.)T 1 // (2/0 fo'iyi) , i A" 6/0 

2 (9/o(2/:) L/o(2/i) J 9 / 0 (2/0 / 0 (2/i) ‘18 XoO/i) 


1 fi(yi) fa" (Vi) i 1 
1 S/o( 2 /i) /o( 2 /i) ‘ 24 




U 1 I 


2*1 


d 0 - 


/= 


,27' F/ /2V\*“|) 

3 t x T: 2 V(tJ 

i Myi) i /i(yQ / j i 

" 12/0(2/1) 16/0(2/1) v 2 yj 1 

1 rp 2 , 2 2 i T,'_Tr fiv Yi 

40 l_~ 2 __r o 3 2 22 2 \ r, 1 2 V rj J 

1 M (yd 1 1 f 2 ( 3 / 1 ) fo' (yi) | 1 rPj' 1 2 j , 9 j //j_ 

"12/o(2/i) ' 12/o(2/0 /o(2/i) + 40 L 2 1 3 id3 Ti 

T, 


rp n /rp / 

2A 2 M-+2An{^- 


/T/V 

T/" , 

rp t r 
-g 11 

2V_ 4 f2VY1 


Ti 1 

0 Ti 

T: 4( vT 1 y l J 

rp r r 

ll 1 1 

(Ti'Y~ 



T: , I 

\tJ _ 




Uf, Trfjdiyi) /(yd/o/yi) / > 

‘ (V ! ^JDo( 2 /i) /o(yi) /o( 2/0 Ji 


16 


D _if 4 ^ ?v _ 9 ( 2 VY 

P2 ~sL 2 Ti 1 T: Ti ‘ \ T } J 


Po= 


V 1 


2 o= 

2 i = 


: 10 L 2 ' 


- 2 A- 2 U, r '' T ' 


Xf )'3 


2 V" 
T, 


2t % 

rp r rp n /rp /\ 2 

, -2^f~2A 2 ^-P2,l 2 (^) - 

I'm, ir2V_, /?VY1 

7 ! 1 5 T: r, 4 Vyi/_ 


‘ 10 L 2 1 3^ s 
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„ A MVi) i Ml 0 / a 

is Uyi) 20 Myi ) 1 5 

i J) (y~-l VMM 


_2Y\_ 

Tj 

1\’ T r " 


\ A o J 1 | n / x 1 

50 U Vl ) L 2 1 3^ 3 2 ^ 2 T x Ti 2 

1 < 


2y\ 3 - 

2i 


ISO ( 
I 7 / 


Mv i ) r 

/o(2/i) L 

a -f(l a+2 ,,) + A + a A [ 2 (^)-f]- 

-KfXfMf)! 


(tt) (1 &+^)+f'+ 


36 


KK? 0 '-fR 

X¥)(R(f)l 


AJ L 

Ti + 


A **- w 


Wi 


■Ait' — aR +i A^Ai 
At/' —Aic A ~ A 2 A k — A*/ Ah 


■djX — A f:+ i 


(2£kg.Q) 
(2£ifc£5) 
(2g£^4) 

"'-A 2 A k "'-ZA 2 'A k l> -ZA 2 "A k '-A t '"A k 


A h '"=A t+x "~A t A ii "-2A t , A t , -A t "A )t 


AA^A li+1 "-A 2 A.J''-4A 2 'A k "'-6A 2 "A k "~ 

4A 2 '"A k '-AAA k 


(k—2, 3) 


(k—2) 

(2 gk S6) 

(0 A m <£'4) 


By — A 2 
By =2 A 3 A z 

By"=2(A,') 2 +2A 3 A 3 " 

B%= A$A$ 

B % — A 3 ' Ai + A 3 A± 

b 9 =a 4 2 

Bio= A 3 A S 

Finally, 

R. P. fli(c)= R. (c) +F? 2 i2 Cc]-A2„ 2 /2 22 (c) 
(4) Evaluation of M 3 (c): 
p 2 (w— c) 2 


n T 
=M 31 (c)-M 0 2 M 32 (c) 




where 

and 


/ N (ie— c) 2 f 5 ' 2 T fv (w—c) 2 , 

(C)- Jn T dy Jv M-cy dy J Ul T dy 


u 


■ w -X 


(w—c) 2 


s' 2 r» fw-c) 2 


dy\~ y >- w - T ^ d y d y 

m 1 J v J vi - 1 


(a) The integral M 32 (c) is evaluated in much the same way 
as Hn(c) ; that is, 

f*b o rbo f'ljB 

6 3 M 32 (c)= Wb 2 (B)b I Tw B 2 d VB d VB - 

JO Jnn Jo 

a 'bo Mbo y*riB 

VJB~d'f\B I I Tv)b dr} B dr} B -\- 

o J VB J o 

J ^bo rbo f*t)B \ 

^ w B djjB | J Tw B 2 di] B dr\ B \-\- 

a 'ba rbo rnB 

dy B \ Tw B 2 drj B di\ B + 

0 J V B J 1} 

r&e rbo 

I W B 2 dy] B I Tr} B dr) B ~h 

Jo J V B 

rbo rto r^B \ 

4 \ w B di) B \ Tw b drj B dt) B )— 

J o Jt} B Jo / 

/ r b * r b o 
2 c 3 ( w B di) B Tt) B dri B A 

\J 0 Jtjd 

rbo rb* fas \ 

Jo ^ B j Jo dyBj-\- 

rbo rbo 

A dltB Tr\B dl}B 
Jo Jvb 

where b Q has the same meaning as in the evaluation of // 22 (c) 
and where 

T= T(w b ) = T x — [(Tj— 1) -2=1 M 0 2 ] w B -^~ i/ 0 V 


The integrals in il/ 32 (e) are evaluated by numerical integra- 
tion using Simpson’s rule. Values of w B are taken from the 
table in appendix B. 

(b) For convenience, the integral M zi (c) is transformed as 
follows: ; 


where 


M n (c) =M 31 i (c) +M ilx (e) - M ih (c ) 

m “- (c> =/: —y i * j" sDp *j: --f ■ ’- 1 ^ 
"‘‘W-r (ilr *£ * 


;« (w—c) 

(w— c) 2 j T 

>y s T 

It is recognized that 

^Z/2 (on — / 

[W -Y L dy=H 1 (c) 
(w—c) 2 


r 

Jti 


Tliexofora . 






By additional transformations, the. following equations are 
obtained: 

M u (c)=H 1 (c) P(c) — Q(c) 
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where 


n , . f* (w— e) 2 J C T J C v - V" 7 

Q(C)= L -^r-*L&=ip iv J, 


- (v>— c) 2 


T 


or 


f 6 (w— c) 2 J f’ T J f b (w—c) 2 j 

Q( ) 6 3 Jo. 46 T dn J 0 .i» (w-c) 2 4 J, 

The integral Q(c) is eyaluated by numerical integration 
using Simpson’s rule; P(c) is evaluated in the calculation 
of H 21 {c). 

The integral Mu (c) is obtained in exactly the same way 
as Ku(c) and H in (c ) ; that is, 


K - (-s+ «■ l- 2 ^ f)+ 

Finally, 

R. P. M,(c) =R. P. il/ 3!l (c)+ITi(c)[fi 2 12 (c) -P(c)]+ 
Q(c)-_Wil/ 32 (c) 

(5) Evaluation of A r 3 (c): 

r * r t! P , f * , , (w-c) 2 , 

14 J yi fo-io — r 1 

Mo dy J* P rfy J, {w-cY dy ~ 


-Vo 4 dy 
J *1 J ffl 


If — c') 2 

jP 1 (i/ 2 -y) dy 


It can be shown that the second and third integrals are 
identical; therefore, 


where 


iV s (c) — A r 3 i(c) — 2Mfh u (/) -f-A/o 4 A r 3 i(e) 


(w— c ) 2 


/yi y.w— w Jfi J- Jv {vt—c) 

N„M - (5V5P */,' <»■-»> ^ *> 

A Wc>= f* f’^ 

J l/i J yi -*■ 

(a) The integral A^fe) is evaluated by numerical integra- 
tion in a manner similar to PAfc), i? 22 (c), and l/ 32 (c) ; that is, 

A r 3t (c)= f dy P ^ c ) dy P dy 

JFi Jift -P J# 

= ^(Jo 4 fp^’ 7 f 

c 2 f*d, Pp d v fdn) 

Jo Jo Jn / 


j Afj4 (c) — ^3 T 'IPbVjJ.bJ^ P dl7s 

r 6 * rfe 

2c P dij B ic B drjs P dijs-)- 

JO Jo Jtjb 

X b f* f} b rk> \ 

T dr\ B drisj T d-q B j 


where 


T= Ti — [(Pi — 1) il/o 2 ] 


U'B — -pj- VoW 


The integrals in N 2i (c) are evaluated by numerical integration 
in a manner similar to that used in the evaluation of d/ 32 (c), 
and so forth. Most of the integrals will already have been 
evaluated in the calculation of i?i(c), iJ 22 (c), and (c). 
(b) For convenience, N 32 (.c) is broken down as follows: 

to-vW 

d «£ ^^ T £! -’ <v>~v>*v+ 
^ iy to-vW 

Let 


where 


N a (e) — A 7 3! } i (c)-i-A r 32 2 (c) — i\ r 3 2 3 (c) 

JY 32 i (c)=£‘ —^— 2 dy£ (yt _ y)iy 

A W c) = £ ^- 2 dyj; (V ^ {y 2 -y)dy 
A^ S ( C ) = J" dyj '* (y 2 — y)dy 

Now, 

A r 32 s (c)= 

J Si 


>! T 
’ si (w— c; 

Since y 2 — y l =1.0, and 


Plfi f'l/J c) 2 

2 dyj* — jr— l(y 2 —yi) — (y-y L )]dy 


-Df \ C H P J f 11 - (w—c) 2 J 

P(C)= J, (SVW-^J, ~r iy 


it is found that 


AWc)=P(e)-P L (c) 


where 


Vl (w—c) 2 , , 

— (y— 2/0 dy 


and 


fsi T C‘ 

p ' (c) -l (S=JF rf!, X 

1 f b T 

= Pjo,46 (W-C) 2 6<l(!j;C)4 

<?i(^;c)= J’ 6 Y V dv-2c£ P 77 d7,+c 2 1 6 


Pi(c) is evaluated by numerical integration using Simpson’s 
rule. Define 

A r 3 2 2 (c)=J^ ^ 2 ^- [(y 2 — y0 — ^/— 2/01% 
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Cvt T 

Li (vj~c) 2 d y~ Ky - 


(0 


X 


v> ( w—c ) 2 J TT r s 

— Y — dy=H l (c) 


N 


t {c)=Ka(p) [i?i(c)- | — (y— yi)#J 


' ~ W T ^ p(w—c) 2 d v J g f o p(w 2 —2cw+c 2 )dr, Jj dr, 

CVl (nn /A2 1 / rbs pVB Ct)B Cbb ('r}K Cl) a 

I — -jr— (j/— yi)dy=ff (j Q Wb 2 ^bJ o T dr, B — 2c | WbC^bJ, T(ii; B +c 2 J o T drn 


f x Vi 2 

The integral — rrn- (y—yi) is evaluated by numerical integration in exactly the same way as N Si (c). 

Jv i 1 

The integral 1Y 32 (c) is* transformed as follows: 


But 

and 
so that 
where 


yi-Vc= (ys—yi) — (yc—yi) =i—<? 
A 7 32 1 (c) = (1— <r)H 2il {c) —J^ic) 

{y ~ v,) dy 


The mtegral J 2 (c) is evaluated in the same way as K n (c). 
Thus 


p j M .1 fo- 1 , J 4 \ , /d-Afl/.W , */ 4 /,(y,) , 1 yi(ViW(Vi) ) Y1 , 1 ,. r 

• R ^u(o)- a _ (r L4+. J+ln j [4 m) - +* \ 5 *>+4 { ^ TOT J Aj + 12 a + * 

a 5 D 0 -\-cr — 4a 3 Co+« 4 (5Z? 0 +0'i)J+ <P ^+6a 2 (7o— 10a s Z? 0 — 4a s (7 1 ^ + o- 3 4aC 0 +10a 2 i?o+6a 2 (7i^+ 

" 4 [“1 2>o-rf i[g-4af7 1 -5aZ) 0 ]+^ (-§ 2o ) 


where 


Finally 


o -2_o_IiM 

0 5 16/o(yi) 

p _Pt L MXM MvOfoivi) \ 

° l ~5 16 (/o(yi) [/ 0 (y,)] 2 f 

n —So i 4 j 1 Uy.) 

^- 6 + 25 / 0 ( 2 / i ) ^ 20 / 0 ( 2/0 
R. P. PV32 (c) = R. P. tVggj (c) T A i 32 2 (c) — A r 32 3 (c) 
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(c) Define 


a *'< c >= jr (sSy d »si *jr iy 


After several transformations, the integral A 7 3 i(c) is brought into a more convenient form 

(c) = A r 31l (c> + K n ( c ) [2 H 2l j (c) + S 2 ]2 (c) ] — A r 3l3 (c) 


where 


„ , , f* T , f» (w-c) 2 7 f« T 

A + (c)_ J, t (w-c) 2 (fy J Sl T ^J, fxr-c) 2 ^ 


A r su(c) = 


Jr tW-C)- Jr, T \) m 


(w—c) 


2 dy 


The integral iV 3 i (e) is evaluated by numerical integration using Simpson’s rule. Some of the integrations have already 
been performed. The integral is given as 


As ' 3 ^ (w— c) 


r ,i,r 

■j JQAb -i J0.4* V.W 


T 


-c) 


2 


The integral A+^c) is evaluated in exactly the same way as K n (c ) ; that is, 


1 \ 


^ l CO= 7 ^j- 3 /o( yi ) + - 


‘ jjPo/o(l/i) -|/o'(3/i) + + [(Pi-2o)/o(y0 +Po/o'G/i) ~\Sz"iy i)_ j -[In(a-cr) ] 2 5- 


2 In (a — <r) 
sa— <r) 


-a 3 [§2>o/o6/0 +^/i(y>)]+^+<r j^/o'&i) -f /oG/i) + a 2 [?o/o(l/i) +|/i(l/i)_ - 


\E-a*F + + 

O 1 


■|/o(yi)+^yo'(yi)-a |_|/o , (l/i)+ 2 >o/o(yi)-r|/i(yi)J-ra 2 -E+^-i f ’| + 

+ [i?o/o(jfi) - |/o'(lfi) -§/i(tfi) -|/o"(2/0 -a(E+ £)]+* 4 [(Pi-2o)/o(2/i) +Pa/o'(2/i) )+In(a-<r) [•5Ino-+ 

+ ( ~^/i(2/i) +§/."(tfi) -f jfi'foi) + <**£-« { 2>o/i(2/i) +| } +§ ^)+ 

+ [po/.&O -§//&.) [}/."G /0 -PafSto) + (a-P:)/i(2/i)] ]-|/o(yi) - 


9 Uo(yi)+<rfa'(yi)l+ofi(yi)'ki «-*r 

(a — a) 2 cr F 1 


, (g — <6 
6 


^Myi)—¥o(.yi)Po - 


< a ~ c'AroMyj +f^/ (i/Opo ■ 


2 LACyOl 2 


L 


. 3 6 jiv»i^u 27 /oG/0 ' 12-' 

1 _2 1 ^ 

j}Ji(yi)—4pJo(yd— 4 3 j o/o / ( yOj— g (a-cO/iG/i) — g (ffl— cr) 2 In (a — ff)/ t (yj — q (a-a) 2 <r/i'(l/i)la(a-o-) 

(-^/oGfcW |/ 0 (yi) 2 >o+|yi(tfi)]- 


/o(yi)2of + 


( a — O- ) 3 In (a — 0 -) { ^ i/i ( 3 / 1 ) £0 


9 / 06 / 1 ) 

t/o^CJ/i) ’Tgfi'tyi) ~ J [Po/i , (l/i)+:Pi/o(l/i)] + 4 <Zo/oG/i) + 4/2Q/1) ~ i 


4 

lata — tr) 


}> 


— J/iCyO— 


(T 3 


f/FF/i) , 1 L/i(y .)] 2 
J 2 + 9 / 0 (yi) 




4 (|/i"(2/i)+| ~[j^) + I 4/1^1)20+4 2 / o yj 

^=^/(yi) ++ [|//(ffi) —2>o/i (1/.)]++ [|/i"(yi)-Po/i , 0/i) + (2o-p0/i(yi)] 


1 /2(1/11/1(1/1) _ 2 ? [/1 ( 3 /i)l 3 

4 /o(i/i) /j 


where 
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n IMVi)] 2 , _ / ( ..s l/2(yO/i(37i) 

c ~ Po “a mu ,) ' 


D=Pofi'(yi)- 





E=Pofa'(yi) +pJ 0 (yO +| //( 2 / 1 ) 


F=p<>.fi iyi ) + ffo/o ( 2 / 1 ) - 1 / 2 (y0 


EVALUATION OF /if”) 

The functions /* < " 1> ( 2 / 1 ), which appear repeatedly in the 
evaluation of the integrals iv/c), i? 2 (c), and so forth, are 
evaluated in terms of Wi {k) and T x ik) as follows: 


fkVc- 


1 

(w c ')\k-l)k\ 



fo im) (yi) = -(Tg 0 ) v ™ 

= ~[r, ;m) go+ mT^-» gi + m 1 [m ~ 1) T^-»g 2 + ...+ 
__i_ 2 1 I lm_f )<7 r -f , . , (0gmg6) 

where 


1 

50 («?,')* 


g 1 =—2g <] Ai 


92 — — 2 (<7o4l 2 ' + 9iA 2 ) 


g m =-2[g 0 A^ m -^+(ni-l)g 1 A^ m ~ 2) + . ■ .+ 

(m _ r ) ! r l 9rA^Z r -V _|_ - . . + p OT _ 1 J. 2 ] 

( 1 \ On) -r 

=(T'g a ) v ™- f-i foT),™ 

(0 ^ m ^ 5) 


/ 4 «(y0 = ^ (y iv <7 0 +2T'"f 7 O t . 1 Cffl, +j2 

jg (F'g 0 S 3 ) h im) -j2 [Myi)Si] v ™ (0^m^2) 

fvO =4§o (^Vo+I (m> +4 ( y,/ 

^(T''9o-5a),/ ra) +^(TM),/ M) ~[/o(yi)S s ]v/^ 

(O^mgl) 

/ 6 ( yO= § 4o (gig jv"S,+ 

1 TW/Cfi 1 m / c* 1 T Q \ 

^‘+600 ^ ' b5- 3600 il ^V 

where 

ft <*> = m <ft> - 3ft < »> - \ <« (0^//<3) 

ft W) = -| aL 5 <*> + - 1 2C , a (t) + y A <S) + 2ft «> (0 ^ k <> 2) 

ft ® = -| ^ «> + 3ft «> + 5ft w > -O 5 C 4 w - 20 C # <*> - 

ftr~ — y A y\-3B«-\- 6fto — 18Cs+^ ft — GO Og-f - 7 5Z/ — 

f C 7 +150D 6 -225E 3 + 3l ~F s 


where 

| (ftoftift' 1 ^--^ L/o(»i)ft], 1 OB> 

and 

ft (t, ==§ft W) -§A ( *> (0^4) 

/3 Cra, (l/i)=^ (ft"<7o+§ w+i ( T'g 0 S 2 ) y 

^ L/o(2/i)ft] ! , I 0a) . (0^m^3) 


ft (m) and ft (ffl) are defined as previously. 

C' 4 ( * > =i (ft^+ftW-ft^O 

ft (&> =^ (ft (i >+fta (t, -ft ci+1> ) 
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C 7 =B 7 
Os = -t4g.il3.ri 4 

a 8 ™ 

A (t) =| (2C 2 ' t ’+C i w -C 3 lt+1) ) 

zv* J =| (2C , s (i) + c s (t) — cm 1 ' 1 ) 

A (t) =| (2C S <*>+ dV^-CV^ 1 ’) 

£,<*>=£),«>— \ Z?,< i+1 > 

4 

^(*>=1 (3i> 5 <*> +ZV*> -A (l+1) ) 

0 

ORDER OF MAGNITUDE OF IMAGINARY PARTS OF INTEGRALS Hi, 

Mu AND Nz 

In the detailed stability calculations the contributions of 
the imaginary parts of the integrals H z , M z , N s , and so forth, 
to the function v(c) are considered to he negligible in com- 
parison with the contribution of the imaginary part of 
Ki(c). A calculation of the orders of magnitude of I. P. H 2 (c), 
I. P. M 3 (c), and I. P. A r s (c) from the general expressions 
given in the preceding pages shows that this step is justified, 
at least for the values of phase velocity c that appear in the 
stability calculations. 

For example, 

I.P. HM = T-P- H hi {c)=*A{w c yUy c ) 


where 


A=-~ 


Therefore 


3 T c (w/) 


3 . 0(ri) + 


c 3 


I-P-^C^-^WW) 


The contribution of I.P. H 2 (c) to r(c) is approximately equal 

tcc 

3 T {w ') ’ w ^- ere I- P- Ki(c). The quantity in 

the brackets is of the order of 0.03, at most, in the calculations 
of the present paper. (In the approximate calculations of 
Re for Mach numbers very much greater than unity, c 

becomes large because cj>l— tj; however, a is small when 

Aio 

c is not much greater than 1 — ttt- and the results of the 

XV Jo 

calculations of R a based on the approximation 

CT mtti 

v(c) P. K 7 {c) are qualitatively correct (fig. 7).) 

h 

From the expression for N z (c), 


I.P. iV,(c) 


2{Wt'Y 


I-P. K n (c) 


so that the contribution of I.P. N z (c ) to »(c) is approxi- 
F a 2 c 2 "1 

mately equal to r 0 \ Ww 7 j 2 J' The < l uail ^ty bi brackets is of 

the order of 0.06 at the most. 

The imaginary part oO/ 3 (c) is considerably smaller. In fact, 

I.P. M 3 (c) ~ Q T f (Wi y I-P- KM 

and the contribution of I.P. MM to o(c) is approximately 
r — c 6 a 2 1 

equal to c 0 /p i’ The quantity in brackets is of the 

order of 0.001 at maximum c. 



APPENDIX B 


CALCULATION OF MEAN-VELOCITY AND MEAN-TEMPERATURE DISTRIBUTION ACROSS BOUNDARY LAYER AND THE 

VELOCITY AND TEMPERATURE DERIVATIVES AT THE SOLID SURFACE 


The mean-velocity and mean-temperature profiles for the 
several representative cases of insulated and noninsulated 
surfaces are calculated by a rapid approximate method that 
gives the slope of the velocity profiles at the surface with a 
maximum error of about 4 percent in the extreme case, for 
which T^O.70 and M Q =0.70. The surface values of the 
higher velocity derivatives and the temperature derivatives 
required in the stability calculations are obtained directly 
from the equations of mean motion in terms of the calculated 
value of the slope of the velocity profile. The Prandtl 
number is taken as unity. 

MEAN VELOCITY-TEMPERATURE DISTRIBUTION ACROSS BOUNDARY 

LAYER 


In a seminar held at the California Institute, of Tech- 


nology in 1942, the. present author has shown that a good 
first approximation to the mean velocity distribution across 
the boundary layer is obtained by assuming that the vis- 
cosity varies linearly with the absolute temperature. With 
this assumption, the velocity w({) is the same function of 


the nondimensional stream function £=— ===== as in the 


Blasius case, 
surface i)—y* 

when <r= 1. Actually, the approximation w(t)=vj B (Z) is 
the first stage of an iteration process applied to the differ- 
ential equations of mean motion in the laminar boundary 
layer, in which po:T l ~‘ (e is a small parameter equal to 0.24 
for air), and w(r)=ia B (f)+ evjx(f) + . . . Calcu- 
lation of n\ (f) for Ti— 1.50 and !Ti=2.00 for fi/ 0 — >0 showed 
that the iteration process is rapidly convergent; the con- 
tribution of the second term to the slope of the velocity 
profile at the surface is 5 percent for T 1 = 1.50 and 8 percent 
for I 7 ! =2.00. In the present calculations the maximum 
error in the slope introduced by taking w(f)=Wzs(f) is about 
4 percent in the extreme case. (See reference 15, in which 
the authors make use of a linear viscosity-temperature 
relation. See also reference 23.) 

That for a linear variation of viscosity with 

absolute temperature is seen directly from the equations of 
mean motion in the laminar boundary layer. The equation 
of continuity is automatically satisfied, by taking 


a nd t he corresponding distance from the 
/ JfoT ; i s obtained by a simple quadrature 

f 


and 



Po* 


c )f* 
dx* 


The stream function \p* and the distance along the surface 
x* are selected as independent variables following the pro- 
cedure of Von Mises, and the dynamic equation o_f mean 
motion becomes for zero pressure gradient 


42 5u*_5 

p °' bx*~b^* 


' du*\ 

u* P* «* 


Define the nondimensional stream function f by the relation 


r= 


v A 

form: 


v Q * « 0 * x * 


The dynamic equation takes the following 


t dw d ( 

'2 3F“3Fv 


dw 


Since p= 7 p in the boundry layer, if y= T, the dynamic equa- 
tion in this form is identical with the equation for the 
isothermal Blasius flow, that is w(f) =vj b ({), or the value of 
the velocity ratio w is equal to. the Blasius value at the same 


value of Jr. The corresponding value of v—V* 


Uq * 

V 


, the 


nondimensional “distance” from the surface, is obtained as 
follows: 

M 


or 


pVM = ^ v * “°* x * by 

dt 

pW =T v 


r*- r 

Jo pVJ Jo 


r rpdt 
0 to 


If (r= 1, the energy and dynamic equations have a unique 
integral and 


T=T t - (T t -I)-^Mo 2 


w- 


7—1 


MJxo 1 


as shown by Crocco. Therefore, 

' =Tl Jl iHV'-u-hr «•’] f- 1 ? 1 


But ~w(t) =w s (f), an( l 


p* ___5i {/* 
Po* u *~by* 



w B di) B -- ~~ 
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f ® i B 

The integrals I w B d-q B and w B 2 d-q B are given in the fol- 
lowing table, and the mean-velocity and mean-temperature 
profiles can be calculated rapidly by this method. (The 

values of are used in the approximate calculation of 

Re . (appendix C).) 

cr ml a 


Vb 

WB 


J>**« 

/Ski's 
\ dg / b 

0.00 

0.0000 

0.0000 

0.0QQ0 

0.3320 

.20 

.0064 

.0066 

.0003 

.3319 

.40 

.1328 

.0265 

.0024 

.3314 

.GO 

.1989 

.0599 

.0081 

.3300 

.SO 

.2647 

.1065 

.0189 

.3274 

1.00 

.3298 

.1660 

.0367 

.3230 

1.20 

.393S 

.2385 

.0630 

.3165 

1.40 

.4563 

.3236 

.0993 

.3079 

1.60 

.5168 

.4210 

.1468 

.2967 

1.80 

.5748 

.5302 

.2064 

.2825 

2.00 

.6298 

.6508 

.2792 

.2663 

2.20 

.6813 

.7821 

.3654 

.2483 

2.40 

.7290 

.9231 

.4643 

.2280 

2.60 

.7725 

1.0733 

.5776 

.2064 

2. SO 

.8115 

1.2319 

.7034 

.1835 

3.00 

.8-460 

1.3978 

.8411 ' 

.1618 

3.20 

.8761 

1.5702 

.9897 

. I40S 

3.40 

.9018 

L74S0 

1. 1478 

.1180 

3. GO 

.9233 

1.9306 

1.3145 

.0936 

3. SO 

.9411 

2.1171 

1.4SS4 

.0805 

4.00 

.9555 

2.3067 

1.6682 

.0640 

4.40 

.9759 

2.6933 

2. 0419 


4. SO 

.9878 

3.0863 

2.4280 


5.20 

.9942 

3.4S2S 

2. 8211 


5.60 

.9975 

3. 8812 

3.2LSQ 


6. 00 

.9990 

4.2805 

3.6167 



With the approximation that y varies linearly with the 
absolute temperature, the slope of the velocity profile at the 
solid surface is simply related to the slope of the Blasius 
profile. Thus 


since w(t) =w B (£), 


and 


or 


c )w dw d'$ 
dtj ~ d£ dij~ pW 


dw 

d? 



dw 

dy~ 


, 0.332 , 

■■Wi — —ji — o 


than the Blasius value and only 2 percent lower at a Mach 
number of 2.0 (fZ\= 1.S1). For the noninsulated surface, 
with 71=0.25, the value of the skin-friction coefficient at _ 
A/ 0 =0 is only 7 percent greater than the Blasius value and 
12 percent, greater at a Mach number of 3.00. 

Since the shear stress at the surface is unchanged in first 
approximation, the boundary-layer momentum thickness has 
the same value as for the Blasius flow 

=0.6667 

V r 0 *x* 

The expression for the displacement thickness 5* gives a 
measure of the effect of the thermal conditions at the solid 
surface and the free-stream Mach number on the thickness 
of the boundary layer. By definition, 

5 = f (1 - P w)d n 
V v a *x* do 

From the relation between dr; and drj B 

pUT-D-Hl -ws)]d ns 

V v • Jo 

= 1.73 + (T 1 -l)1.73-f^A7 0 2 (0.6667) 

= 1.737 T l + I =^ A/ 0 2 (0.6667) 

For the Blasius flow 



1.730 


The “thickness” of the boundary layer b is given by 
6 = 5.60 + ( T 7 ! - 1) 1 . 73 + (o .6667) 


and the form parameter is 


where b is the value of y at the “edge” of the boundary layer 
(when w reaches an arbitrarily prescribed value close to 
unity). It is seen that the shear stress at the surface (or 
the skin friction) has the same value as in the Blasius case 


/ 5 «*\ — . — - /bw\ £>jj 

^ ) = 11, * Un* I — I ^~ L 




V>y 


Pl* Wo* 




= W)*Mi W 


by 


Pi 




The reliability of this approximation can be judged from the 
calculations of the skin-friction coefficient in reference 24, 
in which ^ccT’ 0 - 79 . From figure 2 of reference 24, the value 
of the skin-friction coefficient for an insulated surface at a 
Mach number of 3.0 (7 T 1 = 2.823) is only 12 percent lower 


fi = 2.50T I + 


-Mi? 


For the insulated surface, 

77=2.50+3.50 

CALCULATION OF MEAN- VELOCITY AND MEAN-TEMPERATURE 
DERIVATIVES 

Because of the sensitivity of the stability characteristics 
of the laminar boundary layer to the behavior of the quantity 

dy t ^' e va ^ ues ^he re< I u ^ :re< ^ velocity and tempera- 

ture derivatives at the surface are calculated directly from 
the equations of mean motion, with y=T’ x (m= 0.76 for air). 
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Now ^=f'=5pw so that the dynamic equation is 

—b~w' = {T m v}')'. Since f(0)=f / (0)=0 > 

m 


TiW 

-t i 


where T x 


/ =Wi / j^-^ 2 — M 0 2 — (Ti— 1)J, if cr= 1 . In other words, 

the value of Wi" is readily computed from the value of w/. 
In general, w^ k) is determined from the relation 


— ^ (i;w') k -*=(T m w') k - 1 


or 


(*) 


Wi w = 


„mr/ „_ n , (&—!)! (T m )i" ,, n , 

+ ( F-3)i2! T t » ^ a ' 2> +--- + 

(t-l)l (T% 


0» 


(k— 1— s)s! ^ 
b 




(y^jW-n 

T r 


Wi - 


_b_ ~ 
2T]“ _ 


(/c— 2}f 1 < t-3) 'UJj/ / + . . . -f 


(A"- 2)! 


{k— 2 


^l!_. f(t _ 2 _ r)Wi(r+1)+i + (A- 3)^-2j f//wi( ,_ 3) - 


where 


= 6 (p— l)p l ( - I, ~ 2) w l '+. . .+ 


— -^Ezill „<?-«- 

(p— i— i)!g! p 


and 


P = 1, 2, ... 5 


Pr 


'T 


T,' 


Pi y 2 


Pi 


2(T/) 2 T/ 


T, 3 


rp /. 

„ >n « V / l 

Pi "til ^->3 


„ (TV) 3 ZV" 

■b ^2 


(T“) 


1) yr+m ^ 




(T m ) / 


(2Y) 


Td 


(m, 2) ±-7pz-+3m(m— 1) 


T/Ti 


y 2 “f“ ^ 


2*, 


iI31A { „ i3) ffiT +6(mi2) <3V)UV 


T,” 


Ti 3 


^2^[4T 1 , Td / '+3(T 1 '0 2 j+m ^ 


(Z!I v =(m 4) iZVI 5 +10fm 3) 


T, 4 r 

[15(T, / )(Ti'0 2 +10(T 1 0 2 T/"] + 

m(m — 1) 


(T“) v 


Ti' 


Ti 2 

= (m, 5) 


10Ti"T, / "+5T 1 'Ti W ]+m 


Ti T 


(TiO 


T, f 


■15(m, 4) 


(TiTTf 


[45(T 1 ')*(T l ") a +20(T 1 ')*T 1 '"] + 
[60T ! / T 1 "T 1 /// +15(T/) 2 Ti lv +15(T I / ') a ] + 
- ( -- y7 1} [10(T I ' , 0*H- 1 5 Ti" T^+ G T/Ti T ] + m T f 


where 


(m, 2) =m(m— 1) (m— 2) 
(m, n) = m(m— l)(m— 2) 
m = 0.76 

(m, 1) = — 0.1824- 
(m, 2) = 0.226176 
(m, 3) = -0.506634 
(m, 4) = 1.641495 
(m, 5) = -6.959939 
T/—aw i 


(m—n) 


a=^^ 0 2 -(Ti-l) 


T/^awr-ii-VMJiw/y 

T l '"=aw L " , ~3(y-l)M 0 2 w l 'w l " 

Ti w =aWi v —(y—l)Mo 2 [3(Wi ") 2 +4w I / w 1 "'] 

Ti v = aw ?— 5 (7 — 1 ) M? (2to 1 / 'w/ "d- Wi / «> 1 tv ) 

Td 1 = c Wl VI - ( T - 1) il/ 0 2 [l 0 (Wi " 2 + 1 5Wi "Wi lv + 6w/ w, v ] 

Each velocity derivative is determined from the knowledge 
of all the preceding derivatives. 



APPENDIX C 


RAPID APPROXIMATION TO THE FUNCTION (l-2X)c(c) AND THE MINIMUM CRITICAL REYNOLDS NUMBER 


In section 5, a criterion was derived for the dependence of 
the minimum critical Reynolds number R e on the local 

, . ~ CT min 

distribution of mean velocity and mean temperature across 
the boundary layer. It was found that 

R 6 [T(Cq)] 1 - 76 

>CT min T x C;?- VI — lW(l — Co ) 2 

where c 0 is the value of c for which (I — 2X)o(c) = 0.580 and 



A rapid method for the calculation of the function (1— 2\)v(c) 
and the minimum critical Reynolds number is developed by 
making use of the approximation that the viscosity varies 
linearly with the absolute temperature (appendix B). 
(Since the effect of variable viscosity on the mean-velocity 
profile is overestimated in this approximation, the values of 
Rt (fig. 6(a)) calculated bv this method are lower than 

the values calculated for [s=T°- 7i when heat is added to the 
fluid through the solid surface and higher when heat is 
withdrawn from the fluid.) 

For n=T, the dynamic equation (appendix B) is 

t bw 
~2 bit 

and therefore 

T : b (\ bw\ 
bw by \TbiiJ 
brj 



But 

bT^l fbT\ 
bn T\bn ) B 

so that 

T 2 b f 1 bw C . 2 fbT\ 

bwbn\T bn / j_ 2 ' T\bij ) B _ 

brj 



where 


T=T l - [(T t -1) -^1 



Finally, 



The required values of w B , 


fbw\ 


and £ are obtained from 


the table in appendix B. 

The small correction to the slope X(c) is easily calculated 
once the mean velocity profilehasbeen obtained (appendixB). 
Thus 


X(c)=^P|-l 0.(0) =0) 


The quantity (1 — 2X)r(c) has been calculated as a function 
of c for various values of T x at M 0 = 0, 0.70, 1.30, 1.50, 2.00, 
3.00, and 5.00, and the results of these calculations are 
given in the following table. The decisive stabilizing 
influence of withdrawing heat from the fluid at supersonic 
velocities is illustrated in figure 7. 


! » 

Cl 

Rt„ 

Ti . 

Ct 

Rt , 

cr wt>* 

II 

o 

3/* =0.70 

0.70 

0.1945 

3650 

0.70 

0. 1670 

8440 

.80 

.2695 

1QS0 

.80 

.2390 

2110 

.90 

.3485 

402 

.90 

.3265 

613 

1.25 

.5435 

67 

1.25 

.5425 

74 

1.50 

.6240 

36 

1.50 

.6265 

3S 

Afa=1.30; c>0.231 

c>0.333 

0.90 

0-2455 

9230 

1.30 

G.3450 

2770 

1.05 

.4075 

392 

1.35 

.4585 

275 

1. 20 

.5170 

121 

1.40 

.5505 

99 

1.3422 

.5450 

92 

1.4556 

.6276 

49 

1.50 

.6355 

42 

L60 

.7732 

16 

if»=2.00; c>0.500 

i/«=3.00; e>0.667 

1.63 

0. 5074 

671 

2-4S 

0.6730 

186 

1.65 

-543S 

207 

2.52 

.7058 

59 

1.70 

.6155 

75 

2. 62 

.7655 

24 

1. 75 

.6749 

40 

2.72 

.8105 

14 

1.81 

.7275 

25 

2.77 

.8295 

10 

1.85 

.7612 

19 

2.8225 

.8500 

9 

3ie=o.OO; c>O.SOO f 




5.19 

a goes 

174 




5.20 

.8036 

80 




5.30 

.8262 

23 




5.75 

.9008 

6 




6.G625 

.9350 

3 
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APPENDIX D 


BEHAVIOR OF 


FROM EQUATIONS OF MEAN MOTION 


In order to determine the effect_of free-stream Mach 
number, thermal conditions at the solid surface, or 
free-stream pressure gradient on laminar stability, 
it is necessary to know the relation between . these 
physical parameters and the distribution of the quantity 


p ~j— across the boundary layer. The value of 


d ( dw\ 

TyVTy) M 


solid surface is obtained directly from the dynamic equation 

/ , • a \ i mi i . p d~ f (AaX , , i 


(equations (74) and (79)). The value of 


/ dw\ 
x dy) 


at the 


surface, which is also useful in the discussion of laminar sta- 
bility, is obtained from the dynamic and energy equations 
as follows : 


~d 2 ( dw\l rd 2 /wy 

jy 2 v p <WJrw wl 


Vh'" 2w 1 ,, T I / wSTS' , 2w l '(T 1 , y 
~~TT -—fr- tt—- 

Differentiating the dynamic equation once yields the result 


/// 2mTiV)i' ... , T 

! - — ^ L 


('V A 2 rp ft- 

m(m— 1) 


At the solid surface the rate of change of temperature 
dT* 

and the rate at which the work is done by pressure 

gradient u* both vanish, and the rate at which a fluid 

element loses heat by conduction equals the rate at which 
mechanical energy is transformed into heat by viscous 
dissipation. The energy equation becomes 

At* bW\] /bttfV 2 
J>y*\ dy*) J, Ml \by*)i 


T 1 " = - cr (y-l)M o 2 (w 1 y-m i -^<0 
Utilizing the expression for W\" and T\" gives 

[*(?)] 


(Wi 

<r(l+m)(y-l)ilf 0 2 ^f 


where 


YA. ( w 'W r > L_ i? 1 

[dy\T)]r ts iiWi 2y+»;^r 


From this, expression for ^ 1C following con- 

clusions, which are utilized in the stability analysis, arc 
reached; 

When|~^(j7^J vanishes, the quantity is still 

positive. 

When the free-stream velocity is uniform, 

&(?)l =<T(1+7n)(T_1)M ° 2 Wl ' 

that is, (yr) * s alwa ys positive. 

When the surface is insulated, 

[ j ? ( t )], - ' k1 +“)(•'- w ~ t )‘ 

and [!•(?)! is always positive, regardless of the pressure 
gradient. 


372 



APPENDIX E 

CALCULATION OF CRITICAL MACH NUMBER FOR STABILIZATION OF LAMINAR BOUNDARY LAYER 


For thermal equilibrium the rate of heat conduction from 
the gas to the solid surface balances the rate at which heat 
is radiated from the surface. If the rate at which heat is 
withdrawn from the fluid reaches or exceeds a certain critical 
value at a given local supersonic Mach number, the laminar 
boundary-layer flow is stable at all Reynolds numbers. 
(See section 6B.) The purpose of the following brief calcu- 
lation is to determine the equilibrium surface temperatures at 
several Mach numbers and compare these temperatures with 
the critical temperatures for laminar stability. (See fig. 8.) 

When the solid surface is in thermal equilibrium, 

f> (|p)d*= [ 0 £ d?[OT*)4- (2?)«]<fe (El) 


where e is the emissivity, <r is the Boltzmann constant, and 
the other symbols have already been defined. (See refer- 
ences 14 and 15.) Consider the case in which the free stream 
is uniform and the temperature is constant along the surface. 
For (7=1, 



r P ^ 



y l 

where stagnation temperature T s equals 

Also ( j = 8 if the approximation i±= T is 

employed. (See appendix B.) Since ki*=c p n 0 *ii ] ,=c p (! Q *Ti. 



= 0.332 c, -yV Po* flo* T 0 * (T s — TO yhWo 



When the integrations in equation (El) are carried out, the 
following relation is obtained for the determination of the 
equilibrium surface temperature: 


where 


V^(TA-1) = (T I -T 1 )VA4 

_ 

C? Po* Po*-\t(y—l)c p T 0 * 


The equihbrium surface temperature under free-flight 
conditions is affected principally by the variation in density 
Po* with altitude h. The results of calculations carried out 


for altitudes of 50,000 and 100,000 feet are given in the 
following table: 





T*—T ier 

(ft) 

Aft 


(flg. 8)~ 

5GXIG 3 

3.G 

0.370 

Q. 355 

IQ0XL0 3 

2.0 

.220 

. 1S5 


In this table, 

rp rp Q 

ls ll cr / fW \ 

(c f ^)ko* T,*^R L 

where 

Q heat withdrawn from fluid per second per unit width of 
surface 

C f skin-friction coefficient for one side of surface 
E «o *L 

r 0 * 

L length of surface 
2\* free-stream temperature 

k a * heat-conduction coefficient of gas at free-stream tem- 
perature 

In these calculations the following data are used: 

«*=0.50 
L= 2 ft 

2M=400° F abs. 

<r=4.SOX10~ 13 Btu/sec/ft 2 /(°F abs.) 4 
c f —7.7‘3 Btu/slug/°F abs. 

^*=■3.02 X10 -7 slugs/ft-sec 
o 0 *=9S0 ft/sec 

7?=3.61X10- 4 slugs/ft 5 at 50,000 ft 
=3.31X10-* slugs/ft* at 100,000 ft 
K= 3.35 X10“ 4 at 50,000 ft 
=3.66X10" S at 100,000 ft 

Since T,— — T lcr for M a ~ 3 at 50,000 feet 

altitude and for M 0 ~2 at 100,000 feet altitude, the laminar 
boundary layer is completely stable under these conditions. 

It should be noted that under wind-tunnel-test conditions 
in which the model is stationary, these radiation-conduction 
effects are absent, not only because of reradiation from 
the walls of the wind tunnel but also because the surface 
temperatures are low — generally of the order of room 
temperature. 
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TABLE L — AUXILIARY FUNCTIONS FOR CALCULATING THE STABILITY 
OF THE LAMINAR BOUNDARY LAYER FOR INSULATED SURFACE 


c 

X 

V 


Hi 

n i 

M, 

N't 

3/0=0 

0.0372 

0.0000 

0.0004 

0.0102 

0. 5220 

0.28S9 

0.C6S9 

0.2399 

.0744. 

.0001 

.0029 

.02S5 

.4748 

.2710 

.0604 

.3004 

.1115 

.0603 

.0099 

.0561 

.4303 

.2590 

.0530 

.3124 

.1486 

.0006 

.0235 

.0010 

.3887 

.2433 

.0460 
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.1857 

.0012 

.0162 

.1430 

.3499 

.2278 

.0103 

.3211 

.2226 

.0021 

.0802 

.2040 

.3139 

.2120 

.0350 

.3230 

. 2594_ 

.0033 

.1284 

.2782 

. 2SGS 

.1958 

.0301 

.3217 

.2060 

.0050 

.1937 

.3670 

.2505 

.1797 

.0256 

.3174 

.3323 

.0071 

.2794 

.4721 

.2232 

.1639 

.0217 

-30S4 

.3682 

.0098 

.3896 

.5960 

.1987 

.1487 

.0180 

.2935 

.4037 

.0131 

.5286 

.7418 

.1770 

. 1350 

.0139 

.27GS 

.4143 

.0142 

‘ 

. 5767 

.7904 

.1711 

.1312 

.0125 

.2618 

A/o=0.50 

0.0362 

-0.0000 

-0.0004 

-0. 0148 

0. 5122 
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0.0443 

0. 1927 

.0723 

-.0000 
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-.0234 
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-.0244 
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. 2280 
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.0007 ■ 
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.2166. 
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.1662 
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.2420 

.2525 
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.0627 

.2807 

.1530 

.0217 

. 2425 

.2882 

.0036 

.1342 

.1103 

.2513 

.1390 

.0138 

• 240G 

.3237 

.0054 

.2010 

.1695 

.2246 

.1247 

.0158 

.2333 

.3588 

.0076 

.2882 

.2412 
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.1104 

.0128 
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.3936 

.0103 
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.1790 

.0963 
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TABLE I. — AUXILIARY FUNCTIONS FOE CALCULATING THE STABILITY 
OF THE LAMINAR BOUNDARY LAYER FOR INSULATED SURFACE— 
Concluded 


1 c 

X 

V 

L 

Hi 

a 

Mi 

at 

— - • 

3 /#= 0.90 

0 . 0334 

Q . 0000 

- 0. 0015 

- 0.0503 

0 . 4S16 

Q. 1303 

0.0180 

0. 0908 

.0667 

—.0001 

—.0047 

-.0972 

.4422 

.1298 

.0185 

' .1133 

. 1001 

- 0002 

0GS2 

-.1339 

.4049 

.1281 

.0185 

.1366 

.1335 

-.0002 

-.0102 

-. 1746 

.3696 

.1253 

.0182 

.1594 

. 16 m 

-.0001 

-.0090 

-.2034 

.3365 

.1213 

.0175 

.1825 

.2002 

.0001 

-.0029 

—.2250 

.3054 

.1163 

.0166 

.2055 

.2335 

.0006 

.0098 

—.2387 

.2765 

. 1103 

.0157 

.2252 

.2666 

.0012 

.0312 

-.2441 

.2497 

.2251 

.1030 

.0143 

.2439 

.2997 

.0022 

.0634 

-.2407 

.0947 

.0128 

.2597 

.3326 

.0034 

.1086 

-.2281 

.2026 

.0855 

.0110 

.2703 

.3652 

.0051 

.1697 

-.2063 

. 1S23 

. 0759 

.0090 

.2674 

.3976 

.0072 

.2496 

—.1730 

.1641 

.0656 

.0060 

.2515 

.4296 

.0098 

.3518 

-.1302 

.1480 

.0560 

.0021 

. 21S5 

.4612 

.0130 

.4805 

—.0784 

.1340 

.0464 

-.0036 

.1431 

.4636 

.0132 

.4913 

-.0744 

. 1330 

.0463 

—.0040 

.1373 

.4512 

.0153 

.5788 

—.0421 

.1261 

.0418 

—.0075 

.1004 

3 /*= 1.10 

0.0990 

- 0.0003 

- 0.0140 

- 0 . 2037 

0.4026 

0. 0673 

0. 0012 

0. 0806 

. 1320 

-.0004 

-.0206 

-.2630 

.3682 

.0686 

.0038 

- 106S 

. 1650 

-. 0005 

-.0255 

-.3166 

. 335S 

.0683 

.0051 

.1319 

. 19S0 

-.0004 

-.0272 

—.3640 

.3054 

.0667 

.0058 

. 15SS 

.2309 
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-.0232 

-. 40*19 

.2770 
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.1864 

. 263 S 

.0002 

-.0125 
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.0581 

.0062 

.2101 

.2965 

.0009 

.0072 

-.4680 

.2263 

.0516 

.0058 

.2293 

.3292 

.0018 

.0382 
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• .2040 

.0431 

.0047 

.2416 

.3616 

.0031 

.0829 

—.5086 

.1837 

.0333 

.0031 

.2454 

. 393 S 

.0049 

.1442 

5239 

.1655 

.0213 

.0005 

.2310 

.4246 

.0097 

.2247 

— . 5516 

. 149S 

.0081 

- 0032 

.1834 

.4572 

.0098 

.3300 

-.5675 

.1350 

-.0060 

— . 0CS7 

.0764 

.4536 

.0126 

. 44Q7 

-. 6112 

.1245 

-.0203 

— . 01o / 

-.0737 

.5104 

.0160 

.5789 

-.6875 

.1151 

-.0360 

-.0230 

-.2366 

M , = 1.30 

0 . 2541 


- 0 . Q561 

- 0. 5982 


0.0244 

0.0003 

0.2200 

.2858 

™BVrrrM 

-.0505 

-.6508 

.2255 


.0016 

.2440 

.3173 


-.0364 

-. 69S7 

.2041 

.0183 

.0014 

. 26-14 

. 34SS 

nHnTgjfl 

-. 0117 

-.7430 

.1845 

.0109 

.0003 

.2742 


.0021 

- 025S 

-. 7856 

.1667 


-.0016 

.2700 

.4111 


.0790 


.1507 


-. 004S 

. 22S5 

.4415 

.0057 

. 1503 

— . SS34 

.1366 

-.0236 

-.0090 

.1184 

.4721 


.2449 


.1242 

— . 0404 

-.0169 

-.0818 

.5020 

.0114 

.3652 

- 1 . 0977 

.1136 

-- 0C2S 

-.0294 

-. 4943 


.0120 

.3893 

— L 1334 

.1119 

-.0671 

-. 0324 

-.5971 

.5416 

.0167 

.5777 

- 1 . 3074 

.1020 

-.0834 

-.0549 

-i. soso 


TABLE IL— AUXILIARY FUNCTIONS FOR CALCULATING THE STABILITY 
OF THE LAMINAR BOUNDARY LAYER FOR NONINSULATED SURFACE 


C 

X 

s 

L 

Hi 

a 

A/s 

AT 

. Vs = 0 . 70 ; 77 = 0.70 

0.0262 

0. 0056 

0.0825 

0. 0635 

0.6102 

0.3272 

0.0524 

Q . 274S 

.0521 

.01 12 

.1645 

.0949 

.5725 

.3157 

.0502 

.2920 

.0777 

.0166 

.2466 

.1184 

.5367 

.3645 

. Q4SI 

.3081 

.1030 

.0220 

.3297 

.1400 

.5026 

.2936 

. 045S 

.3233 

.1281 

.0274 

.4146 

.1632 

.4703 

.2828 

.0433 

.3380 

.1529 

.0327 

.5023 

.1904 

.4396 

.2724 

.0412 

.3519 

.1701 

.0365 

.5661 

.2130 

.4191 

.2651 

.0396 

.3610 

.1726 

.0370 

.5754 

.2163 

.4162 

.2642 

.0394 

.3623 

- 114 = 0 . 70 ; 77 = 0.80 

0.0237 

0.0033 

0 . O4S6 

0.0279 

0.5954 

0.2811 

0.0493 

0. 1369 

.0472 

.0066 

.0965 

.0374 

.5620 

.2737 

.0475 

.1504 

.0705 

.0099 

. 1413 

.0430 

.5300 

.2663 

.0457 

.1635 

.0937 

.0132 

.1925 

.0482 

.4994 

.2590 

.0437 

.1763 

. 1168 

.0164 

.2417 

. 0550 

.4701 

.2514 

.0417 

.1882 

.1397 

.0197 

.2926 

.0649 

.4420 

.2439 

.0397 

.2001 

.1625 

.0230 

.3457 

.0789 

.4152 

.2363 

.0378 

.2110 

. 1S51 

.0263 

.4017 

.0982 

.3897 

.2287 

.0359 

.2213 

.2075 

.0297 

.4614 

.1236 

.3654 

.2210 

.0339 

.2311 

.2298 

.0331 

.5253 

.1562 

.3424 

.2133 

.0321 

.2400 

.2409 

.0349 

.5592 

. 1754 

.3313 

.2094 

.0310 

.2443 

.2475 

.0359 

.5801 

.1877 

.3248 

.2071 

.0303 

.2465 




£ 

II 

© 

o 

17 = 0.90 




0.0433 

0.0036 

0.0517 

0. 0051 

0.5506 

0.2410 

0.0435 

a 1426 

.0863 

.0072 

.1028 

-.0047 

.4939 

.2304 

. 0404 

.1633 

.1291 

.0108 

. 156S 

—.0111 

.4414 

.2191 

.0370 

.1846 

.1714 

.0145 

.2173 

-.0079 

.3930 

.2074 

.0337 

.2032 

.2135 

.0185 

.2885 

.0096 

.3485 

.1951 

.0304 

.2203 

.2551 

.0227 

.3746 

.0462 

.3080 

.1825 

.0272 

. 2339 

.2963 

.0274 

.4805 

. 1071 

.2715 

.1693 

.0240 

.2462 

.3166 

.0299 

.5426 

.1489 

.2547 

.1637 

.0224 

.2517 

.3268 

.0312 

. 5762 

.1776 

.2466 

.1606 

.0217 

.2541 

3/6=0.70; Ti=lJ25 

0.0346 

- 0.0016 

HSBll 

- 0.0476 




0. 1462 

.0692 

-.0032 

-.0476 

-.0797 

.4678 

Knyl 


.1634 

.1010 

-.0048 

— . 069S 

-. 1013 

.4276 

.1661 

.0292 

.1794 

. 13S9 

-.0062 


—.1132 

- 3S3G 



.1956 

.1738 

— . 0076 

-. 1021 

1155 

. 353S 

.1529 


.2108 

.2088 

— . 00S7 

-.1085 

-.1081 

» I ■ 

.1448 



.2439 

-.0095 

— . 10oi 

-.0912 

.2888 

.1354 


.2342 

.2789 

0101 

-.0917 

-.0645 

.2597 

.1249 

.0135 

Koni 

.3138 


-.0641 

-.0281 

■1 ...il 

. 1133 

.0161 

.2409 

.3485 


-.0203 

.0179 

. 20S6 


.0139 

.2297 

.3831 

-.0092 

.0427 

.0734 

.1865 

. 0S70 

.0113 

.2069 

.4174 

-.0079 

. 1286 

.1373 

.1668 

.0728 


.1616 

.4512 

-.0059 

.2414 

.2071 

.1495 

.0582 

.0042 

. 081 $ 

.4846 

-.0031 

.3859 

.2770 

.1345 

.0427 

-.0012 

-.0601 

.5092 


.5184 

.3212 

.1248 

.0314 

-.0067 

2262 

.5190 

.0006 

,o//9 

.3349 

.1212 

.0269 

-.0091 

-.3028 
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TABLE III— PHASE VELOCITY, WAVE NUMBER, AND REYNOLDS NUMBER FOE NEUTRAL SUBSONIC DISTURBANCE (STABILITY LIMITS) 

FOR INSULATED SURFACE 


c 

or 

i? 

as 

Ri 

C 

a 

R 

as 

Rt 



Mo=0 




-_r 

.Mo =0.90 



0. 0372 

0. 0321 

25,500,000" 

0.0038 

3,030,000 

0.0334 

a. 0107 

111,000,000 

0. 0012 

12, 600, 000 

.0744 

.0685 

1,500,000 

.0082. 

178,000 

.0667 

.0248 

5, 960,000 

.0028 

679, 000 

.1X15 

. 1103 

278,000 

.0131 

33,100 

.1001 

.0421 

1, 030, 000 

.0048 

117, 000 

.1486 

.1585 

83,000 

.0189 

9, S80 

.1335 

.0632 

290,000 

.0072 

33,000 

.1857 

.2146 

32, 600 

‘ '.0255 

3,880 

.1669 

."0885 

106,000 

.0101 

12,100 

.2226 

.2808 

14, 800 

.0334 

1,760 

.2002 

.1186 

.46,000 

.0135 

5, 240 

.2594 

.3590 

7, 700 

.0427 

917 

.2335 

IJL540 

22, 600 

.0175 

2,570 

.2960 

.4535 

4,420 

.0540 

526 

.2666 

\ 1961 

12, 100 

.0223 

1,380 

.3323 

.5707 

2,760 

.0679 

329 

.2997 

.2459 

7,020 

.0280 

800 

.3682 

.7243 

1,850 

.0862 

220 

.3326 

.3053 

4,320 

.0348 

492 

.4037 

.9589 

1,360 

. 1142 

162 

.3652 

-.3777 

2, 820 

.0430 

321 

.4143 

1.0770 

1,280 

<1282 

153 

.3976 

,4638 

1,950 

.0529 

222 

.4143 

L2730 

1,530 

-.1515. 

182 

.4296 

15733' 

1,410 

.0653 

161 

.4037 

1.2940 

' 1,880 

.1540 

223 

.4612 

.7260 

1,090 

.0827 

125 

.3682 

1.1960 

3,530 

.1424 

421 

.4636 

.7396 

1,0S0 

.0843 

123 

. 3323 

1. 0400 

-6. 710 

.1238 

.799 

. 4812 

.8810 

1,010 

. 1004 

115 

.2960 

.8728 

13, 300 

.1039 

1, 580 

.4812 

1.0130 

1,230 

.1154 

140 

.2594 

.7177 

27, 500 

■ 0S54 

3,270 

.4636 

1.0120 

1.740 

.1153 

199 






.4612 

1.0070 

1,820 

.1148 

207 






.4296 

.9027 

3, 180 

. 1029 

363 



Mo=0.50 



.3976 

..1823 

5, 590 

. QS92 

637 






.3652 

.6642 

9, 940 

.0757 

1,130 






.3326 

7.5523 

18,500 

.0629 

2, 100 

0.0362 

0.0251 

36,6(50,000 

0.0029 

4,270,000 






.072$ 

•. 0538 

2, 130, 000 

.0063 

248,000 






.1085 

.0868 

392,000 

.0101 

45, 700 



AA=1.10 



.1446 

.1250 

116, 000 

. 0146 

13, 500 






.1806 

.1695 

44, 500 

.0198 

5,190 






.2166 

.2216 

20, 200 

.0258 

2,360 

0.0990 

G.00S6 

5, 730, 000 

0. 0009 

618,000 

.2525 

.2-829 

10,400 

.0330 

1,210 

.1320 

.0268 

769,000 

.0029 

82,900 

.2882 

.3556 

5,850 

.0414 

682 

.1650 

.0468 

224,000 

.0050 

24, 100 

.3237 

.4442 

3, 570 

.0518 

416 

.1980 

.0707 

85,000 

.0076 

9, 160 

.3588 

.5549 

2,330 

.0647 

272 

.2306 

.0991 

38, 300 

.0107 

4, 130 

.3936 

.6993 

1, 620 

.0815 

189 

‘ .2638 

.1329 

39,300 

.0143 

2.030 

.4280 

.9301 

1, 230 

.1084 

144 

.2965 

.T727 

10,600 

.0186 

1,140 

.4306 

.9558 

.1,220 

.1114 

142 

.3292 

.2200 

6,260 

-.0237 

675 

.4362 

1.0140 

1,190 

.1182 

139 

.3616 

.2755 

3,920 

.0297 

423 

.4362 

1.1880 

L 410 

.1384 

164 

.3938 

.3417 

2,610 

.0308 

281 

.4306 

1.2150 

1,580 

.1416 

184 

.4246 

.4159 

1.850 

.0448 

199 

.4280 

1.2150 

1,660 

' .1416 

.. 194 

.4572 

.5193 

1,350 

. 0560 

146 

.3936 

1.1240 

3,080 

.1310 

359 

.4836 

,.6268 

1,100 

.0676 

119 

.3588, . 

.9788 

5, 670 

.1141 

661 

.5104 

.8010 

991 

.0864 

107 

.3237 

.8272 

10,800 

.0964 

1.260 

.5104 

.9165 

1, 220 

.09SS 

131 

.2882 

.6869 

21, 100 

.0800 

2,460 

.4836 

.8927 

2,060 

.0962 

223 






.4572 

.8023 

3,320 

.0865 

358 






. 4246 

. 6785 

5,930 

.0732 

C39 



Mo=0.70 



.3938 

".5706 

10, 400 

.0622 

1,120 

0. 0353 

0. 0191 

53, 40t), 000 

0.0022' 

6, 100, 000 


- - 

Mo= 1.30 



.0705 

.0415 

3,060,000 

.0047 

349. 000 






.1058 

.0677 

555, 000 

.0077 

63,400 






.1410 

.0984 

161, 000 

.0112 

IS, 400 

0.2541 

0.0451 

63,800 

0.0047 

6. 630 

.1762 

.1344 

61, 100 

.0154 

6,9S0 

.2858 

.0818 

24, 800 

.Q0S5 

2,370 

.2114 

.1766 

27,300 

.0202 

3,120 

.3173 

.1202 

12,300 

.0125 

1,280 

.2464 

.2268 

13, SCO 

.0259 

1.580 

.3488 

\ 1636 

6,990 

.0170 

726 

.2813 

.2857 

.7, 630 

.0326 

872 

.3800 

.2132 

4, 280 

.0222 

445 

.3181 

.3570 

4, 550 

.0408 

520 

.4111 

.2707 

2,800 

.0281 

291 

.3505 

• . 4433 

2,900 

'.0506 

331 

.4418 

.3377 

1,930 

.0351 

201 

.3847 

.5515 

; 1,960 

.0630 

224 

.4721 

.4100 

1,420 

.0133 

147 

.4185 

.6951 

.1,420 

.0794 

162 

.5020 

_5123 

1,110 

.0532 

115 

.4352 

.7917 

’ 1,230 

.0904 

141 

.5072 

<5316 

1,070 

.0552 

m 

.4452 

.8655 

. 1,160 

.0989 

132 , 

.5416 

.7582 

886 

. 078S 

92 

.4559 

.9704 

' '1,110 

.1108 

127 

.5416 

.8931 

. 1,080 

.0928 

112 

.4559 

1.1230 

' 1.330 

.1283 

152 

.5072 

.7781 

2,310 

. 0S09 

241 

.4452 

1.1420 

1, 650 

.1304 

189 

.5020 

.7592 

2,550 

.0789 

265 

.4352 

1.1230 

1,980 

.1283 

227 

.4721 

..6458 

4,500 

.0671 

468 

.4185 

1. 0720 

2,670 

.1225 

305 

.4418' 

:5401 

7, 9S0 

.0561 

829 

.3847 

.9381 

4,810 

.1072 

550 






.3505 

. 7965 

8, 8S0 

.0910 

1,010 






.3161 

.6659 

16, 700 

.0761 

1,910 
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TABLE IV.— PHASE VELOCITY, WAVE NUMBER, AND REYNOLDS NUMBER FOR NEUTRAL SUBSONIC DISTURBANCE (STABILITY LIMITS! 

FOR NONINSULATED SURFACES 


c 

a 

R 

ai 

Ei 

c 

a 

R 

at 

Ei 

Mi =0.70; ri=0.70 

3/t=0.70; 77=0.90 

0.0262 
.0521 
.0777 
.1030 
. 1281 
.1529 
.1701 
.1726 
.1726 
.1701 
.1529 
.1281 
.1030 
.0777 
.0521 

0.0339 
.0734 
.1188 
. 170S 
.2308 
.3030 
.3670 
.3 m 
.4986 
.4977 
.4732 
.4175 
.3460 
.2620 
.m3 

82. 400. 000 

5.360.000 

1.110.000 

371.000 

161.000 

83.400 
57,200 

54.400 
69,000 
73,900 

121,000 

270.000 

711.000 
2,500,000 

14. 600.000 

0.0041 

.ooss 

.0143 

.0205 

.0277 

.0364 

.0441 

.0454 

.0599 

.Q59S 

.0568 

.0502 

.0416 

.0315 

.0206 

9.900.000 

644.000 

133.000 
44,600 
19,300 
10,000 

6.870 
6,540 
S.2SQ 

8.870 
14,500 

32.400 

85.400 

300.000 

1.750.000 

0.0433 
.0863 
.129L 
.1714 
.2135 
.2551 
.2963 
.3166 
.3268 
.3268 
.3166 
.2963 
.2551 
.2135 
.1714 • 
.1291 
.0863 

0.036S 
.0815 
. 1353 
.1996 
.2775 
.3728 
.49S0 
.5814 
.6347 
.7817 
.7701 
. 7307 
.6275 
.5133 
.3972 
.2858 
.1793 

17, 100,- 000 
1, 040, QCQ 
200,000 

62.500 

25.500 
12,400 

6,970 

5,520 

4,990 

6,500 

7,920 

11,600 

25,200 

60,300 

170.000 

617.000 
3,740,000 

0. 0042 
.0092 
.0153 
.0226 
.0314 
.0422 
.0563 
. G65S 
.0718 
.0SS4 
.0871 
.0827 
.0710 
.0581 
.0449 
.0323 
.0203 

1,930,000 
118,000 
22.700 
7, 070 
2,880 
1,410 
7S9 
624 
565 
735 
895 
1,310 
2,850 
6,820 
19,200 
-69,800 
423,000 

Mj=0.70; Ti=0.S0 

3/a =0.70; 27=1.25 

0.0237 

.0472 

.0705 

.0937 

.1168 

.1397 

.1625 

.1851 

.2075 

.2298 

.2409 

.2475 

.2475 

.2409 

.2298 

.2075 

.1851 

.1625 

.1397 

.1168 

.0937 

.0705 

0.0237 
.0504 
.0804 
. 1138 
.1509 
.1923 
.2382 
.2908 
.3510 
.4237 
.4668 
.4962 
.630S 
.6233 
.6056 
.5609 
.5062 
.4465 
.3S27 
.3164 
.24S9 
.1822 

157,000,000 

9.910.000 

1.970.000 

633.000 

263.000 

129.000 

70.900 
42,600 

27.500 
18,800 

15.900 

14.500 

18.500 

21.400 
27,200 

44.900 

77. 400 

141.000 

280.000 
630,000 

1. 690.000 

5.890.000 

0.002S 

.0059 

.0094 

.0133 

.0176 

.0224 

.0278 

.0339 

.0409 

.0494 

.0544 

►057S 

.0735 

.0726 

.0706 

.0654 

.0590 

.0520 

.0446 

.0369 

.0290 

.0212 

18,300,000 

1,150,000 

230.000 
73,700 

30.600 
15,000 

8,260 

4,960 

3,200 

2,190 

1,860 

1,690 

2,160 

2,500 

3,170 

5,230 

9,010 

16.400 

32.600 

73.400 

197.000 

686.000 

0. 0346 
.0692 
.1040 
.1389 
.173S 
. 2QSS 
.2439 
.2789 
.3138 
.3485 
.3831 
.4174 
.4512 
.4346 
.5092 
.5190 
.5190 
.5092 
.4846 
.4512 
.4174 

0.0160 
.0346 
.0564 
.0819 
.1120 
.1477 
. 1S99 
.2403 
.3002 
.3722 
.4594 
.5668 
.7061 
.9067 
1. 1800 
L44S0 
1.5880 
L7250 
1. 5370 
1.2580 
1. 0330 

78, 800,000 
4,380,000 

770.000 

217.000 
78,900 
34,000 
16,500 

8,830 
5,070 
3, 110 
2,020 
1,380 
1,000 
760 
643 
615 
640 
806 
1,390 
2,740 
5,360 

0.0016 
.0036 
.0058 
.0QS4 
.0115 
.0152 
.0195 
.0247 
.0308 
.0382 
.0471 
. 05S2 
.0725 
.0931 
.1211 
.14S6 
.1630 
.1770 
.1577 
.1291 
.1060 

8,090,000 
450,000 
79,000 
22,200 
3,100 
3, 490 
1,700 
907 
520 
319 
207 
142 
103 
78 
66 
63 
66 
83 
142 
281 
550 




